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To compare the metal-insulator transitions (MITs) of V0 2 and V2O3 we analyze the relations between the 
structural and electronic properties of the vanadium Magneli phases. These materials set up the homologous 
series V„02 n -i (3 < n < 9) and have crystal structures comprising typical dioxide-like and sesquioxide- 
like regions. As the MITs of the vanadium Magneli phases are accompanied by structural transformations, 
we are able to discuss the effects of characteristic changes in the local atomic environments. The systematic 
investigation of the transport properties is based on a new and unifying description of the crystal structures 
of the Magneli phases including VO2 and V2O3. Our results lead to a comprehensive understanding of the 
MITs in the Magneli class and shed new light on the role of particular electronic states for the MIT of V2O3. 
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1 Introduction 

The vanadium oxides [1-3] comprise compounds with different formal vanadium valency stages, which 
reach from two in VO, three in V2O3, and four in VO2 to five in V2O5. In addition to these configurations 
mixed valent compounds can be synthesized. Amongst the latter the so-called Magneli phases, defined by 
the general stoichiometric formula 

V„0 2n _i = V 2 3 + (n - 2)V0 2 where 3<n<9, (1) 
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Compound V„C<2n-i 
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Table 1 Formal V 3d charges as well as transition temperatures in the series V„02n-i • 



are of special interest since they give rise to a homologous series of compounds with closely related crystal 
structures. This type of homologous series has been reported for the first time by Magneli for the molybde- 
num oxides [4]. Nowadays additional Magneli series are known for the vanadium, titanium, niobium, and 
tungsten oxides. The first structural characterization of vanadium Magneli compounds by means of x-ray 
investigations traces back to Andersson and Jahnberg [5] in 1963. As equation Q indicates, the Magneli 
phases take an intermediary position between V2O3 and VO2, thus between the valency stages three and 
four. In addition to the chemical relationsship, the crystal structures of the Magneli phases actually consist 
of rutile and corundum-type blocks and consequently show structural affinity to both the dioxide and the 
sesquioxide. In section[5]we study the structural relations between these materials in detail - discovering 
the potential of gradually transferring the crystal structure of VO2 into that of V2O3 by making use of the 
Magneli phases. This observation allows for a comprehensive understanding of all these compounds. 

As a function of temperature each Magneli phase undergoes a metal-insulator transition (MIT), except 
for V7O13, which is metallic at all temperatures. The MITs are of first order and accompanied by distinct 
structural transformations. Tablensummarizes transition temperatures, inferred from electrical resistivity 
measurements by Kachi et al. [6], and formal V 3d charges, varying from two in V2O3 to one in VO2. In 
an ionic picture a metal atom contributes three 3d and two 4s electrons, whereas an oxygen atom accepts 
two 2p electrons. As the Magneli phases are characterized by vanadium atoms in mixed valent states, the 
electronic properties were expected to be influenced by charge ordering. Although this is interesting from 
both the experimental and theoretical point of view, only few studies are reported in the literature. 

FigureniHustrates the behaviour of the MIT temperatures listed in table^as a function of the vanadium- 
oxygen ratio. We recognize a broad minimum centered at the ratio corresponding to V7O13. Interestingly, 
the MITs of the vanadium Magneli phases are coupled to an anomaly in the magnetic susceptibility, which 
closely resembles the characteristics of a transition from a paramagnetic to an antiferromagnetic state [6]. 
Nonetheless, the anomaly is not accompanied by long range magnetic order. Yet, all Magneli phases enter 
an antiferromagnetic groundstate at sufficiently low temperatures; Neel temperatures are given in figure^ 
In some sense the anomaly may be interpreted as a precursor effect. There is a maximum of the transition 
temperature for V2O3 and a gradual decrease on approaching VO2, where antiferromagnetic ordering is no 
longer observed. V7O13 is the only deviation from this general trend. Although a continuous reduction of 
V2C<3-like regions in the crystal structure is a possible explanation for decreasing transition temperatures, 
this kind of mechanism has not been confirmed. While MIT and magnetic ordering are coupled for V2O3, 
they appear at different temperatures and thus are likely to occur independently for other Magneli phases. 
As we are interested in the MITs, the subsequent electronic structure calculations do not take into account 
the low temperature antiferromagnetism. Actually, it would not be possible to include the magnetic order 
since the detailed magnetic structures of the ordered phases are not known. 

Specific heat data for the Magneli phases by Khattak et al. [7] revealed values larger than expected from 
pure lattice contributions. Due to the smallness of the magnetic entropy increase at the Neel temperatures 
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Fig. 1 Comparison of the MIT temperatures 
and the magnetic ordering temperatures of the 
Magneli phases. Besides V7O13 all members 
of the series undergo an MIT. Except for VO2 
each compound develops an antiferromagnetic 
order at low temperatures. In the case of V2O3 
the MIT and the magnetic ordering coincide at 
a temperature of 168 K. 



the authors argued in favour of magnetic correlations above the phase transitions. To explain their findings 
they proposed a model based on linear antiferromagnetic chains. However, low temperature investigations 
of the magnetic susceptibility by Nagata et al. [8] ruled out this model. The authors suggested interpreting 
insulating V§Oi5 m terms of a charge-density wave, as Gossard et al. [9] speculated in the case of V4O7. 
They attributed small magnetic moments at the metal sites of insulating V4O7/V6O11 to singlet pairing. 

Single crystal x-ray data by Marezio et al. [10] indicated chains of V 3+ and V 4+ ions in the insulating 
phase of V4O7, whereas the vanadium valences are disordered in the metallic phase. This kind of evidence 
for charge localization arises from the comparison of experimental V-0 distances in V4O7 with estimated 
values referring to an ionic picture. VC>6 units with large/small V-O bond lengths are interpreted in terms of 
y3+/y4+ con fjg Ura tions, respectively. However, due to the covalent portion of the interatomic bonding the 
calculated charges might be poor. Specific heat data by Griffing et al. [11, 12] implied that metallic V4O7 
and V7O13 are characterized by a coexistence of itinerant and localized V 3d electrons. Nuclear magnetic 
resonance experiments confirmed this result and additionally found two inequivalent metal sites, but were 
in conflict with pure S = 1/2 and S — 1 spin states, corresponding to d 1 and cP configurations [9, 13]. In 
the insulating phase the Magneli compounds show indications of increased charge localization and singlet 
pairing, except for V3O5. The charge localization may split the subbands associated with the inequivalent 
metal sites and hence may pave the way for a Mott-Hubbard transition. 

While the MIT temperature of VO2 increases on the application of pressure, V2O3 reveals the opposite 
trend. The Magneli phases reflect the behavior of the sesquioxide as their transition temperatures decrease 
with growing hydrostatic pressure [2] . Moreover, Canfield et al. [14] reported on unifying aspects emerging 
in the Magneli series on the application of pressure. They observed for both V§Oi5 an d (Vo.98Cro.o2)sOi5 
a transition from a paramagnetic metal (PM phase) to an antiferromagnetic metal (AFM phase) and further 
to a paramagnetic insulator (PI phase). Consequently, the PM-AFM transition present for V7O13 at atmo- 
spheric pressure should not be regarded anomalous for the Magneli series. While the PM-AFM transition 
is found in VgOi5 only above 9kbar, it is always present for the chromium doped system. The magnetic 
transition exhibits the same pressure dependence in all three cases (—0.75 K/kbar). Since the MIT temper- 
ature decreases at a larger rate, the AFM modification is stabilized for sufficiently high pressures. Although 
in other Magneli phases the PM-AFM transition is not observed up to 20 kbar, the authors assume different 
pressure dependences of the MIT and the antiferromagnetic ordering as a general feature, common to each 
Magneli phase. This again indicates that the transitions evolve independently. 
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2 Crystal structure and electronic properties of V0 2 and V 2 3 

Above 340 K VO2 crystallizes in the rutile crystal structure based on a simple tetragonal lattice with space 
group P^/mnm (D\i), as shown in figure|2] Using single crystal experiments McWhan et al. [15] deter- 
mined the tetragonal lattice constants an — 4.5546 A and cr = 2.8514 A and the positional parameters. 
With respect to the vanadium sites the VO2 structure is based on a body-centered tetragonal lattice, where 
each metal site is surrounded by an oxygen octahedron. Octahedra at the center and the corners of a rutile 
unit cell are rotated by 90° around the rutile c-axis. This reduces the lattice symmetry from body-centered 
tetragonal to simple tetragonal and causes the unit cell to contain two formular units. Equally adjusted oc- 
tahedra form chains along the rutile c-axis. Within the chains neighbouring octahedra share edges, whereas 
corner and center chains are connected via corners. Next to the filled oxygen octahedra the VO2 structure 
comprises the same amount of empty octahedra, which likewise form chains along the rutile c-axis. Hence 
VO2 is based on a regular three dimensional network of oxygen octahedra partially filled with vanadium 
atoms. Each octahedron has orthorhombic symmetry, but deviations from tetragonal and even cubic sym- 
metry are small. One observes two different V-O distances, where the apical distance appears twice within 
each octahedron and the equatorial distance is found four times [16]. 

We define local coordinate systems centered at the metal sites. Because of the different orientations of 
the oxygen octahedra two reference systems are required. In both cases the local z-axis is oriented along 
the apical axis of the (local) ocahedron, which is either the (110) or the (110) direction. Compared to the 
traditional alignment of the local x and y-axis parallel to the metal-ligand bonds, these axes are rotated by 
45° around the local z-axis. Thus they are parallel and perpendicular to the rutile c-axis, respectively. This 
definition of a local rotated reference frame is useful not only for rutile VO2 but also for related structures. 
Referring to the original rutile coordinates we introduce the pseudorutile axes a pru t, b piu t and c prut , which 
allows us to discuss the crystal structures of the different Magneli phases on a common basis. 

Figuref3]displays the angular parts of the metal d orbitals relative to the reference frame of the central 
metal atom. The cubic part of the crystal field splitting results in a separation of the d orbitals in threefold 
degenerate ti g (d x 2_ y 2, d xz , d yz ) and twofold degenerate e g {d?, z 2_ r 2, d xy ) states. While both e g orbitals 
point to oxygen atoms, the d x 2_ y 2 orbital points along the rutile c-axis (local x-axis) and the local y-axis. 
Moreover, d xz ld yz orbitals are directed to faces of the local oxygen octahedron. Hence a and 7r-type V-V 
overlap along the rutile c-axis is mediated by d x 2_ y 2 and d xz states, respectively. Due to the 45° rotation 
of our coordinate system the d x 2_ y 2 and d xy orbitals are interchanged compared to the standard notation. 
The lobes of the d yz orbital point perpendicular to the rutile c-axis. Because of the ratio chIcir = 0.6260 
of the crystal axes they exhibit reduced er-type overlap, compared to the d x 2_ y 2 states. While both d x 2_ y 2 
and d yz orbitals connect metal atoms separated by tetragonal lattice vectors, coupling between metal sites 
at the center and the corners of the unit cell is mediated by d xz orbitals, which exhibit remarkable overlap 
with the d x 2_ y 2 states belonging to the vanadium atoms of neighbouring octahedral chains. 

The first order phase transition of stoichiometric VO2 at 340 K is accompanied by structural distortions 
leading to the monoclinic Mi phase, which is based on a simple monoclinic lattice with space group P2\/c 
(C| h ) [17]. Despite very similar gross features, distinct differences distinguish rutile and monoclinic VO2. 
First, a characteristic metal-metal pairing along the c prut -axis is present in the monoclinic structure, which 
modifies the V-V distances (2.851 A in the rutile case) and gives rise to alternating values of 2.619 A and 
3.164 A. Second, zigzag-type displacements of the metal sites evolve along the diagonal of the rutile basal 
plane (local z-axis). The shift direction alternates along a pru t/c pru t but not along 6 prut . Due to the zigzag- 
type distortions two different apical V-O bond lengths of 1.77 A and 2.01 A are observed. In addition, the 
V-V pairing results in two short (1.86 A and 1.89 A) and two long (2.03 A and 2.06 A) equatorial distances 
in each oxygen octahedron. Third, a lattice strain is present in the monoclinic configuration. 

Uniaxial stress along the rutile (110) axis or doping with some percent of chromium/aluminium yields 
two additional insulating phases: the monoclinic M2 and triclinic T phase [18-20]. In the Mi phase V-V 
pairing and zigzag-type lateral displacements of the metal atoms affect each vanadium chain. In contrast, 
in the M2 phase half of the chains dimerize and the other half show zigzag-type deviations. The T phase 
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Fig. 2 Crystal structure of metallic VO2. The in- 
serted cuboid determines the unit cell of the rutile 
structure with the c-axis oriented upwards. Large 
and small spheres represent vanadium and oxygen 
atoms, respectively. The vanadium atoms form in- 
finite chains along the rutile c-axis. Furthermore, 
each vanadium site is surrounded by an oxygen oc- 
tahedron. The latter are oriented uniformly in each 
vanadium chain but rotated by 90 ° with respect to 
the rutile c-axis among neighbouring chains. Hence 
they establish local coordinate systems. 

is intermediate since the dimerized M2 chains start to tilt and the zigzag chains evolve a dimerization until 
eventually the Mi structure is reached. 

An instability of the rutile crystal structure is found not only for VO2 but also applies to other transition 
metal oxides, as for instance M0O2 and NbC>2. For the latter compounds the relations between structural 
distortions and electronic structure changes have been studied in [22,23]. As VO2, they are characterized 
by metal-metal pairing and lateral zigzag-type displacements in the low temperature phase. Despite all the 
similarities a striking difference regarding the in-plane distortions of the metal atoms and the modifications 
of the surrounding oxygen octahedra distinguishes monoclinic VO2 and M0O2. Only for the latter material 
the octahedra follow the in-plane metal shifts, whereas for VO2 the vanadium atoms move relative to their 
oxygen octahedra. This relative shift affects the bonding between V 3c? and O 2p states as it influences the 
orbital overlap. Thus an antiferroelectric distortion of the VOg octahedra is inherent in the monoclinic (Mi) 
phase of VO2. From the different behaviour of VO2 and M0O2 it was concluded that the antiferroelectric 
distortion cannot explain the tendency of the transition metal oxides to form distorted variants of the ideal 
rutile structure. For the same reason the mechanism must be excluded as a possible driving force for MITs 
connected to the destabilization of the rutile crystal structure [24]. 

Theoretical approaches towards understanding the MIT and the simultaneous structural transformations 
range from Peierls [25,26] to Mott-Hubbard [27,28] schemes. Possible driving forces are lattice instabili- 
ties, electron-phonon interaction, and electronic correlations. The importance of the lattice for stabilization 
of the different phases of VO2 and the symmetry change at the transition point to strong electron-phonon 
coupling [15]. This line of reasoning is supported by a comparison of rutile and monoclinic VO2 via ultra- 
sonic microscopy, which reveals a strong elastic anisotropy in the metallic but almost no anisotropy in the 
insulating phase [29]. The electronic structure of metallic VO2 has been probed by optical measurements 
revealing the lowest empty V 3d levels 2.5 eV above the top edge of the O 2p bands [30]. Ultraviolet and 
x-ray photoelectron spectroscopy show a 8.5 eV wide occupied band directly below the Fermi energy [31]. 
The valence band splits up into low and high binding regions with widths of 1.5 eV and 6 eV, respectively. 
While the low binding contributions are attributed to the V 3d states, the broader part of the valence band 
mainly traces back to O 2p states. According to oxygen K-edge x-ray absorption spectroscopy unoccupied 
V 3d bands extend from the Fermi energy to 1.7 eV and from 2.2eV to 5.2eV [32]. For monoclinic VO2 
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Fig. 3 Angular parts of the transition metal 
3d states with respect to the local coordinate 
system of the central metal atom in the dis- 
played rutile unit cell. Note that the rutile c- 
axis is oriented upwards. The t^g manifold, 
resulting from the crystal field splitting of the 
metal states, comprises the d x 2_ y 2, d xz and 
d yz orbitals, whereas the e g manifold com- 
prises the d 3z 2_ r 2 and d xy states. While the 
d x 2_ y 2 lobes point at edges of the local octa- 
hedron, the d xz and d yz orbitals are directed 
towards faces. Both e 9 orbitals point at oxy- 
gen atoms. 

photoelectron spectroscopy shows a sharpening accompanied by an energetical downshift of the occupied 
V 3d bands. Shin et al. reported a band gap of 0.7eV for the insulating phase [33]. 

Goodenough proposed an energy band scheme for both metallic and insulating VO2 [1]. Starting with 
electrostatic considerations he placed the O 2p states well below the V 3c? states. Octahedral crystal field 
splitting leads to lower V 3d tig and higher V 3d e g levels, where the former are located in the vicinity of 
the Fermi level and split into du (oriented parallel to the rutile c-axis) and it* states. In the monoclinic case 
metal-metal pairing causes the d\\ band to split into filled bonding and empty antibonding states. Due to 
antiferroelectric zigzag-type displacements of the metal atoms, the 7r* bands shift to higher energies. Shin 
et al. found the du band splitting to amount to 2.5 eV and the ir* bands to rise by 0.5 eV [33]. In contrast 
to the Goodenough scheme Zylbersztejn and Mott proposed an MIT mechanism based on strong electron- 
electron correlations [34]. According to these authors especially the one-dimensional du band is affected 
by correlations rather than by electron-lattice interaction. In the metallic phase correlations are efficiently 
screened by the ir* electrons, but screening is reduced below the phase transition as the ir* bands are subject 
to an energetical upshift. Thus the narrow du bands become susceptible to strong Coulomb repulsion and 
undergo a Mott-transition - opening a correlation gap. Recently, the MIT of VO2 was explained in terms 
of a Peierls instability of the one-dimensional dn band, embedded in a background of the remaining V 3d 
tig states [24, 35]. Molecular dynamics calculations by Wentzcovitch et al. likewise point to predominant 
influence of the lattice degrees of freedom [36]. 

Stoichiometric V2O3 at ambient pressure undergoes an MIT at 168 K accompanied by a transformation 
from the high temperature corundum into the low temperature monoclinic structure, where additionally an 
antiferromagnetic order appears (AFI phase). Systematic studies have provided the temperature, pressure, 
and composition dependence of the MIT in the (y \- T ^x)2^>z system [37, 38]. Doping with chromium or 
aluminum results in a paramagnetic insulating configuration (PI phase), which retains the lattice symmetry 
of the paramagnetic metal (PM phase). Since the metal-insulator boundary terminates at a critical point, a 
gradual PM-PI crossover appears when decreasing the temperature from well above the 500 K region. For 
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Fig. 4 Crystal structure of V2O3. For simplicity the schematic view does not include the oxygen sublattice, which is 
an octahedral network qualitatively equal to the case of VO2. The vanadium atoms once more occupy octahedral sites 
but now form hexagonal structures perpendicular to the Ch cx -axis. Hence vanadium pairs arise both along a pr ut=Chex 
and Cp ru t. The table records experimental V-V distances (in A) for the paramagnetic metallic, the antiferromagnetic 
insulating, and the paramagnetic insulating phase of V2O3 [40,41]. 



a chromium concentration 0.005 < xci < 0.018 three phase transitions are encountered: PM-PI, PI-PM, 
and PM-AFI, where the last two are of first order [39]. For xci < 0.005 the reentrant MIT disappears. In 
contrast, on doping with titanium the PM-AFI transition temperature decreases and for xti > 0.05 only a 
metallic phase is found. Hydrostatic pressure has almost the same effect as doping with titanium. 

Corundum V2O3 is based on a trigonal lattice with space group R3c (D® d ). Oftentimes a non-primitive 
hexagonal unit cell comprising six formula units is used, whereas the primitive trigonal unit cell contains 
two formula units. The metal atoms are octahedrally coordinated by six oxygen atoms, where the oxygen 
network is qualitatively equivalent to the arrangement in VO2. We identify the Ch ox -axis of the hexagonal 
unit cell with the a pmt -axis. Two thirds of the oxygen octahedra are filled with vanadium, but no contin- 
uous chains of face-sharing VO@ octahedra arise. Instead, two filled octahedra are followed by an empty 
one. Neighbouring chains are shifted along Ch ox , giving rise to hexagonal vanadium structures, see figure 
|4] This forms the basis of describing the corundum structure in terms of a hexagonal cell. Because of the 
trigonal symmetry each metal atom exhibits three identical in-plane V-V distances of 2.87 A. Metal-metal 
electrostatic interaction along Ch ox (through octahedral faces) causes V-V anti-dimerization reflected by a 
bond length of 2. 70 A, which is considerably larger than the ideal octahedral separation of 2.33 A. Hence 
the hexagonal vanadium structures deviate from being planar. There are three V-O bond lengths of 1.97 A 
and three other of 2.05 A in each VC<6 unit. 

At 168 K stoichiometric V2O3 distorts into a monoclinic structure with space group 12 /a (CfjJ, char- 
acterized by increasing nearest neighbour V-V distances [40]. Bond lengths across shared octahedral faces 
along a pru t increase from 2.70 A to 2.75 A. Moreover, the threefold degeneracy of the in-plane V-V dis- 
tance is lifted, giving rise to a strongly elongated (2.99 A) bond along c pru t-axis. The symmetry breaking 
at the PM-AFI transition results in two crystallographically inequivalent oxygen sites, where each oxygen 
octahedron comprises four sites of the first and two of the second kind. Although the V-O distance varies 
froml.95Ato2.llA, the average length is similar to the PM phase. In contrast to the symmetry breaking 
at the PM-AFI transition, the corundum structure is conserved at the PM-PI transition [41]. Nevertheless, 
V-V bond lengths across octahedral faces/edges increase to 2.75 A/2.91 A. A comparison of PI (aluminum 
doped) and AFI V2O3, as based on hard and soft x-ray absorption spectroscopy, revealed the same local 
symmetry for the metal sites and hence (on a local scale) identical distortions [42]. 
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V2O3 is generally regarded as canonical Mott-Hubbard system [27,28]. A description in terms of the 
one-band Hubbard model is based on a level scheme by Castellani et al. [43]: Crystal field splitting of the 
V 3d states due to the octahedral coordination yields lower t 2g and higher e g states. Because of the trigonal 
lattice symmetry the former split up into a\ g and e g contributions. The non-degenerate a ig orbital points 
along Chc X and gives rise to a one-dimensional band due to covalent bonding. This electronic level scheme 
was confirmed by band structure calculations [44] . Castellani et al. assumed the bonding a\ g orbital as fully 
occupied, whereas the antibonding states shift energetically above the e g levels. This leaves one electron 
per vanadium site in the twofold degenerate states, making the system susceptible to degeneracy lifting 
distortions. One e g electron (total spin S = 1/2) suggests using the half filled one-band Hubbard model 
as the simplest model describing V2O3 [45]. 

Experimentally, the PM-AFI phase transition was analyzed by photoelectron spectroscopy. As reported 
by Shin et al. [33], the spectra show O 2p states in the energy range from — lOeV to — 4eV, whereas V 3c? 
states occur within 3 eV below the Fermi level. Except for slight modifications of the bandwidth no drastic 
change is found in the V 3d band structure below the MIT. However, the DOS at the Fermi energy is rather 
small even in the metallic phase. A high resolution photoemission study by Shin et al. [46] revealed a band 
gap of 0.2 eV for the insulating phase. The PM-PI transition of (chromium doped) V2O3 was investigated 
by Smith and Henrich [47] using photoelectron spectroscopy. In the insulating phase at room temperature 
they found a low emission intensity at the Fermi level, which increases after cooling into the metallic state. 
High temperature PI spectra reveal a larger V 3d bandwidth than low temperature AFI spectra, which may 
be due to thermal broadening or absence of magnetic order. According to oxygen K-edge x-ray absorption 
spectroscopy unoccupied V 3d bands extend from eV (Fermi energy) to 6 eV with maxima at about 1 eV 
and3eV [48]. 

The model of Castellani et al. has been called into question by the results of polarized x-ray absorption 
spectroscopy, which indicate a vanadium S = 1 spin state [49]. Furthermore, the first excited states miss a 
pure e g symmetry but include a\ g contributions, which requires an explanation beyond the pure one-band 
Hubbard model or beyond models projecting out the a\ g states by means of simple dimerization. LDA+U 
calculations succeeded in explaining the electronic and magnetic properties of the AFI phase, in particular 
the peculiar antiferromagnetic order [50]. While every vanadium spin is parallel to the adjacent spin along 
Chcx, there are two antiparallel spins and one parallel spin at the nearest in-plane vanadium sites. Enforcing 
spin degeneracy in the calculation, influence of the crystal structure on the electronic properties was shown 
to be rather small. Introducing the antiferromagnetic order yields an optical band gap of 0.6 eV. Moreover, 
LDA+U findings point at a spin S = 1 model without orbital degeneracy. In contrast, by means of model 
calculations Mila et al. [51] proposed a S = 1 spin state and orbital degeneracy in the AFI phase. Starting 
with the assumption of strong covalent bonding in the vanadium pairs along Ch cx the authors supposed the 
intersite ai g hopping matrix elements to dominate. However, an analysis of the hopping processes in V2O3 
revealed hopping integrals between second, third, and fourth nearest vanadium neighbours being equally 
important for the shape of the a lg band [52]. Evidence for orbital ordering in the AFI phase is taken from 
magnetic neutron scattering [53]. 

LDA calculations display only minor response of the electronic structure to crystal parameter changes. 
But a slight narrowing of the characteristic a\ g bands in the insulating phase formed the basis for describing 
the PM-PI transition by a combination of LDA calculations and the dynamical mean field theory (DMFT), 
demonstrating the influence of electronic correlations in the insulating phase [54]. Confirming predictions 
of the LDA+DMFT approach, a prominent quasiparticle peak was observed in the photoemission spectrum 
of metallic V2O3 [55]. However, despite many studies the electronic properties and MITs of V2O3 are still 
controversially discussed. The answer of the electronic states to the structural modifications at the MIT of 
VO2 paved the way for interpreting the transition as an embedded Peierls instability. A similar procedure 
seems to fail for the sesquioxide since here the structural modifications hardly influence the LDA results. 
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3 Crystal structures of the Magneli phases: a unified representation 

To analyze the electronic properties of the vanadium Magneli phases, we first have to understand in detail 
the crystal structures - a prerequisite for relating electronic changes at the MITs of the Magneli phases to 
structural transformations occuring simultaneously. We elaborate a new and unifying point of view of the 
atomic arrangements underlying all members of the series. A distinct advantage of this representation is its 
applicability to the crystal structures of vanadium dioxide and sesquioxide as well. Thus a comprehensive 
understanding of the whole Magneli series V n C>2n-i including its end members VO2 (n — ► 00) and V2O3 
(n = 2) is eventually achieved. The detailed modifications of the various crystal structures at the MITs are 
discussed in subsequent sections. 

As mentioned before and as is obvious from the stoichiometric relation 0, the crystal structures of the 
Magneli phases are usually interpreted as rutile-type slabs separated by shear planes with a corundum-type 
atomic arrangement [5]. The rutile-type slabs extend infinitely in two dimensions and have a characteristic 
finite width corresponding to n NOq octahedra in the case of the compound V„02 n -i- Since octahedra at a 
slab-surface share faces with octahedra from a neighbouring slab, the atomic arrangement at the boundary 
is closely related to the corundum structure. The rutile crystal structure is therefore disturbed and adjacent 
slabs are mutually out of phase, thus giving rise to the denotation shear plane. Precise investigation of the 
changes of the local atomic coordination on passing through the Magneli series is rather complicated using 
the ordinary representation of the crystal structures. Thus it is advantageous to describe the structures in a 
different manner. For that purpose we start with the oxygen sublattice, which turns out to be most similar 
for all the compounds. As demonstrated for VO2 and V2O3 the sublattice can be understood as a regular 
space filling network of neighbouring oxygen octahedra. The octahedra are mutually connected via edges 
along Cprut and via faces along a pm t/& pr ut- In order to illustrate this geometrical arrangement schematical 
projections (parallel to a pru t) of sandwich-like O-V-O slabs cut out of the crystal structures of the dioxide 
and the sesquioxide are shown in figure|5] A sandwich-like slab consists of a vanadium layer confined by 
oxygen layers at both ends. Due to the projection the oxygen octahedra appear as hexagonal structures and 
the octahedral network becomes a two-dimensional hexagonal network. 

Both the infinite vanadium chains of VO2 (along c prut ) and the hexagonal arrangement of the vanadium 
atoms in V2O3 are clearly visible in figure|5] Actually, the basic structural features of the dioxide and the 
sesquioxide are captured completely in this graphical representation. Adjacent slabs yield the same pro- 
jection as shown except for a possible translation perpendicular to a pru t. To be more specific, in adjacent 
VO2 slabs filled and empty oxygen octahedra, i.e. filled and empty oxygen hexagons, are interchanged. 
In neighbouring V2O3 slabs the vanadium atoms are shifted by one hexagon along the c pm t-direction. To 
sum up, provided that the oxygen networks of VO2 and V2O3 are exactly the same, the only difference 
between the compounds arises from a different vanadium sublattice. The configuration of the latter is op- 
timally described in terms of the presented slab-type projection. In particular, the infinite vanadium chains 
belonging to the dioxide are contrasted with finite chains of length 2 in the case of the sesquioxide. Such 
a description of the V2O3 structure in terms of 2-chains along c prut instead of hexagonal arrangements is 
preferable for a systematic discussion. There are small differences between the oxygen sublattices of VO2 
and V2O3. However, they do not affect the gross features of the oxygen arrangement, but still are essential 
for a detailed comparison of the materials. As indicated by our representation in terms of O-V-0 slabs we 
find alternating vanadium and oxygen layers along a prut in both compounds. While the oxygen layers of 
the sesquioxide are almost flat, we are confronted with a distinct buckling in the case of the dioxide, which 
is caused by an elongated apical V-O distance. The buckling of the oxygen layers is associated with the fact 
that octahedral faces oriented perpendicular to a pmt in V2O3 adopt a tilt in the VO2 case. Nonetheless, the 
characteristic feature distinguishing vanadium dioxide and sesquioxide certainly is their different pattern 
of filling the oxygen network with vanadium atoms. 

Having developed a common description for the crystal structures of the end members of the Magneli 
series it is easy to extend this point of view to the remaining compounds. For that purpose the structures of 
V4O7 and VgOn are depicted in figure [6] analogous with the representations of VO2 and V2O3 in figure 
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Fig. 5 Projection parallel to a prut of a O-V-O sandwich-like slab cut out of the crystal structures of VO2 (top) and 
V2O3 (bottom). Large and small spheres represent vanadium and oxygen atoms, respectively. Due to the projection 
the oxygen octahedra form a regular hexagonal network. For VO2 infinite chains of vanadium atoms run along the 
pseudorutile c pru t-axis, where the oxygen octahedra are connected by common edges. The projection of an adjacent 
O-V-O slab results in a similar configuration - but empty and filled octahedra are exchanged. Thus the metal atoms in 
VO2 lack nearest neighbours in the a pru t -direction. For V2O3 finite chains of two vanadium atoms run along the c prut - 
axis. The projection of an adjacent O-V-O slab reveals a shift by one octahedral site along c prut . As a consequence, 
the vanadium atoms exhibit exactly one nearest vanadium neighbour in the a prut -direction. 



[5] Except for a slightly different buckling of the oxygen layers, which is visible in figure0for V4O7, the 
aforementioned regular three-dimensional oxygen network forms the basis of all Magneli phases. Hence 
the hexagonal arrangement of the oxygen atoms in the projection of the O-V-O sandwich-like slabs does 
not change. In contrast, each particular material is characterized by the way of filling the oxygen network 
with vanadium atoms. Obviously, there are chains of vanadium atoms running along the c pm t-axis, which 
comprise four and six atoms in the case of V4O7 and VgOn, respectively. In general, there are vanadium 
chains of length n in V„02 n -i- Looking back to the dioxide (n = 00) and the sesquioxide (n — 2) this 
rule coincides with our finding of infinite chains and 2-chains, respectively. In a particular compound all 
metal chains reveal the same length, defined by the parameter n; n — 1 empty oxygen octahedra separate 
them in the c prut -direction. Vanadium chains neighbouring along 6 prut overlap by roughly half the in-chain 
V-V distance, thus giving rise to the characteristic al chain end arrangement depicted in figure[6] This chain 
end configuration is conserved in the whole Magneli series and only the length of the chain center changes. 
For VO2 there are actually no chain ends and the 2-chains of V2O3 contain no center atom. 

The atomic arrangement represented by the projections of the O-V-O slabs is unaltered for slabs neigh- 
bouring in the a prut -direction. Only the vanadium sublattice is shifted along the c prut -axis by an amount of 
n—1 octahedral sites, i.e. n — 1 hexagons in the sandwich-like projection. Thus the last vanadium atom of 
a metal chain is seated on top of the first atom of a chain in the adjacent slab, see figures[8]and[9] Hence the 
chain end atoms exhibit exactly one nearest vanadium neighbour along the a pru t-axis, whereas each chain 
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Fig. 6 Projection along a prut of a O-V-0 sandwich-like slab from the crystal structures of V4O7 (top) and VsOn 
(bottom). Large and small spheres represent vanadium and oxygen atoms, respectively. As in the cases of vanadium 
dioxide and sesquioxide the oxygen octahedra form a regular hexagonal network and the octahedral sites are partially 
occupied by metal atoms, giving rise to vanadium chains along c prut . The chains consist of n — 4 or n — 6 atoms, 
separated by ra — 1 = 3 or n — 1 = 5 empty octahedra. Along & prut the chains overlap by roughly half the intrachain 
V-V distance. The projection of an adjacent O-V-0 slab results in a similar configuration - but the vanadium sublattice 
is shifted along c pru t by n— 1 = 3 or ra— 1 = 5 octahedral sites. As a consequence, chain end atoms exhibit one 
nearest vanadium neighbour in the a prut -direction, whereas center atoms reveal none. 



center atom is surrounded by two empty sites in this direction. All the above considerations apply to each 
Magneli compound. Chain center sites hence mirror the local atomic coordination known from VO2 with 
nearest vanadium neighbours only in the chain direction. Instead, the coordination of the chain end atoms 
resembles the crystal structure of V2O3 with one nearest vanadium neighbour along a prut and c pmt . Due 
to the varying length of the vanadium chains in a qualitatively unchanged oxygen sublattice, the different 
Magneli phases have different stoichiometric V:0 ratios reaching from 1:2 in VO2 to 2:3 in V2O3. As the 
crystal structures of the Magneli phases are connected to each other simply by removing or inserting chain 
center atoms, the Magneli series allows us to transfer a dioxide-type atomic arrangement step by step into 
a sesquioxide-type arrangement. Using these systematics we gain insight into the interplay between VO2 
and V^C^-like regions and their effects on the MIT. Taking into account the mixed valent metal states, the 
Magneli series transfers the electronic configuration of the dioxide {d 1 ) into that of the sesquioxide (d 2 ). 

In the following sections two particular Magneli phases are discussed in detail: V4O7 and VeOxi. To 
define the notation of the different metal sites figure|8]gives a schematic view of the V4O7 crystal structure 
where the regular network of oxygen octahedra has not been included for simplicity. We identify vanadium 
4-chains parallel to the c pru t -direction. Due to mutual interconnection of the chains along a prut a kind of 
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stair-like vanadium arrangement is observable. Stair endings resemble the corundum atomic coordination 
(figure^}, whereas intermediate plateaus are due to the rutile regions of the crystal. The longer the central 
part of the vanadium chains, the larger the separation of the stair endings and the dioxide-like character of 
the crystal. The illustration of VgOn in figure[9]is different from figure[8]due to longer vanadium chains. 
We find four crystallographically inequivalent vanadium sites in the case of V4O7 and six sites for VgOn. 
In both structures these sites set up two types of metal chains: for V4O7 we find the series V1-V3-V3-V1 
and V2-V4-V4-V2, whereas for V 6 O n the series V1-V3-V5-V5-V3-V1 and V2-V4-V6-V6-V4-V2 arise. 
Vanadium layers comprising either chains of the first or of the second kind alternate along a pru t. 

The V-V distances given as insets in the above figures address the dominating structural changes accom- 
panying the MIT of both V4O7 and VgOn [57-60]. At the phase transition of the former compound the 
dimerization in the 1-3 chains becomes stronger and an additional dimerization evolves in the 2-4 chains. 
Thereby the dimerization patterns are reversed because the longer V-V distance appears at the chain ends 
and the chain center, respectively. The V1-V2 bond length is almost constant. For VeOn changes in the 
1-3-5 chains are small although a little dimerization is present at low temperatures. In contrast, the 2-4-6 
chains evolve a pronounced dimerization at the transition. Again the variation of the VI -V2 bond length is 
by far the smallest observed. A profound analysis of the structural changes induced by the MITs is given 
in the subsequent sections where we relate LDA results to local structural properties. However, signatures 
of the V02-like dimerization are reflected by both V4O7 and VgOn. 

As far as the question of choosing the unit cells of the Magneli phases is concerned, several definitions 
have been used in the literature. Andersson and Jahnberg [5] defined primitive translations in terms of the 
rutile lattice (a a, b^, ca) = M(aR, ha, cr), where they distinguished the cases n odd and n even 



-1 1 \ / -1 1 

M odd =11 1 1 , M cvcu = 1 1 1 I . (2) 



re— 3 2 — n 6—n 
2 2 2 



n — 3 2 — 71 6 — 71 



However, this proposal is not useful for a comparison within the Magneli series since for every parameter 
value n the vanadium chains extend in different directions of the unit cell, and due to the two families of 
cell parameters. Le Page and Strobel [61] hence proposed an alternative set of primitive translations with 
the c-axis parallel to the vanadium chains 



h L I = 
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V1-V3 V3-V3 


V2-V4 V4-V4 


V1-V2 


T = 


298 K 


2.97 2.79 


2.93 2.93 


2.77 


T = 


120 K 


3.03 2.67 


2.83 3.03 


2.78 



Fig. 8 Crystal structure of V4O7. For simplicity this schematic view does not include the oxygen sublattice. The 
metal atoms form chains of length n = 4 along Cp rut . There are 4 inequivalent sites in the chains V1-V3-V3-V1 and 
V2-V4-V4-V2. Layers perpendicular to a pru t either comprise chains of the first or of the second kind. The table gives 
measured V-V distances (in A) for both the metallic (298 K) and insulating (120 K) phase of V4O7 [59]. 



While the unit cell of Andersson and Jahnberg is A-centered, the cell of Le Page and Strobel is I-centered. 
Finally, Horiuchi et al. [62] introduced a primitive unit cell with translation vectors 

a H \ ( -1 1 \ / a fl \ 

b ff = 1 1 1 b fl . (4) 

c H J \ n - 5 n - \ I \ c R J 

Because this unit cell has half the size than the previous cells we adopt it for our calculation. Application 
of relation (0} based on the parent rutile lattice is feasible whenever the crystal symmetry is triclinic with 
space group PI {C}) [63], which is the case for the Magneli compounds with parameters n — 4, 9. In 
contrast, V3O5 crystallizes in a body-centered and simple monoclinic lattice with space group 12/ c (C| h ) 
and P2/c (C| h ) above and below the MIT, respectively [64,65]. 

4 Structural and calculational details 

The first vanadium Magneli compound we discuss at length is VgOn, see section |5] In doing so we can 
benefit from a suitable application of the knowledge we already acquired by our former considerations of 
VO2. Consequently, it is convenient to choose VgOn as the starting point of a comprehensive analysis of 
the whole Magneli class. Similar to the MIT of vanadium dioxide and the PM-AFI transition of vanadium 
sesquioxide, the MIT of VgOn at 170 K is accompanied by distinct structural transformation. According 
to an x-ray study by Canfield [60], VgOn crystallizes in a triclinic lattice with space group PI {C}). The 
vanadium and oxygen atoms are located at the Wyckoff positions (2i): ±(x, y, z). For the triclinic lattice 
parameters the author reported a L = 5.449 A, b L = 7.010 A, c L = 31. 437 A, a L = 67.15°, (3 L = 57.45°, 
and 7 L = 108.90° at 298 K. Below the phase transition (at 20 K) he found the parameters a L = 5.495 A, 
b L = 6.944 A, c L = 31.484A, a L = 67.40°, j3 L = 57.13°, and 7L = 108.61°. These values are related 
to the unit cell of Le Page and Strobel as denoted in equation (0, where n = 6. The primitive translations 
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V1-V3 V3-V5 V5-V5 V2-V4 V4-V6 V6-V6 


V1-V2 


T = 
T = 


298 K 
20 K 


2.95 2.91 2.80 2.95 2.85 2.82 
2.90 2.99 2.79 3.06 2.62 3.30 


2.77 
2.77 



Fig. 9 Crystal structure of VeOn. For simplicity this schematic view does not include the oxygen sublattice. The 
metal atoms form chains of length n — 6 along c prut . There are 6 inequivalent sites in the chains V1-V3-V5-V5-V3- 
VI and V2-V4-V6-V6-V4-V2. Layers perpendicular to a pru t either comprise chains of the first or of the second kind. 
The table gives measured V-V distances (in A) for both the metallic (298 K) and insulating (20 K) phase of VeOn [60]. 



of Horiuchi et al. take the form 
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and coordinates given in the L-system can directly be transformed into the i/-system 



L 








V + 




z )l 


V 1/2 



(6) 



The translation vector has to be added because of different choices of the origin in both coordinate systems. 
With the help of equation the positional parameters reported by Canfield for both the high and the low 
temperature configuration of VgOn yield the structural input for the band structure calculation, which is 
summarized in table [2] In addition to the six inequivalent vanadium atoms there are eleven oxygen sites. 
Due to the inversion symmetry of the crystal lattice the unit cell contains 34 atoms. Further structural data 
for metallic VgOn were given by Horiuchi et al. [62]. However, we use the data of Canfield. 

To allow for an adequate interpretation of the electronic structure results, it is important to preserve close 
relation between the representation of the VgOn structure and the parent rutile structure. We arrange the 
unit cell of the L-system in a Cartesian coordinate system such that the alignment of the oxygen octahedra 
resembles the rutile arrangement; we assume sll = (— a\, 0, a%), t>L = (a\, 62, 03), and = (ci, C2, C3). 



An application of the elementary relations a L 



atb 



a 2 T — a\. Furthemore, from the expressions 



b\ + o|, and a L • h L = a L b L cos~/ L 



immediately yields cq, 3 — 1/2 • (a| T- a^L cos7l) and b\ — u, L 
sll ■ cl = a^CL cos/?l as well asb^-CL = &lcl cosai one finds C1/3 = (&lcl cosq;l =F o-lcl cos 0l — 
^2C2)/2a 1 / 3 , which implies a quadratic equation for C2 because c 2 L = c\ + c\ + c|. Solving this equation 
completes the calculation of the unit cell in the L-system and a transformation via equation (|5jl yields the 
unit cell in the H-system. For the high temperature structure of VeOn at 298 K we find (A = 8.6665 as) 



a H = A 



-1.0000 
0.0000 
0.6417 



A 



1.0000 
0.9616 
0.6417 



c H = A 



-0.0053 
5.4629 
3.6352 



(7) 
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High temperature structure 


Low temperature structure 


Atom 


X 


y 


z 


X 


V 


z 


VI 


0.9931 


0.4467 


0.0624 


0.9986 


0.4633 


0.0584 


V2 


0.9731 


0.9264 


0.0672 


0.0081 


0.9568 


0.0612 


V3 


0.9652 


0.4583 


0.2422 


0.9578 


0.4379 


0.2506 


V4 


0.9560 


0.9463 


0.2432 


0.9362 


0.9362 


0.2484 


V5 


0.9765 


0.4689 


0.4170 


0.9885 


0.5086 


0.4074 


V6 


0.9823 


0.9777 


0.4152 


0.0582 


0.0473 


0.4036 


Ol 


0.6888 


0.4114 


0.0332 


0.6584 


0.3744 


0.0402 


02 


0.3202 


0.5570 


0.0666 


0.3173 


0.5650 


0.0658 


03 


0.6845 


0.4670 


0.1090 


0.6809 


0.4685 


0.1098 


04 


0.3736 


0.6180 


0.1460 


0.3464 


0.5937 


0.1508 


05 


0.6930 


0.4060 


0.2134 


0.6954 


0.4057 


0.2150 


06 


0.3138 


0.5883 


0.2468 


0.3360 


0.6176 


0.2396 


07 


0.6697 


0.4485 


0.2968 


0.6683 


0.4469 


0.2974 


08 


0.3171 


0.5339 


0.3474 


0.3161 


0.5435 


0.3468 


09 


0.7039 


0.4223 


0.3898 


0.6826 


0.3998 


0.3940 


O10 


0.3180 


0.5976 


0.4254 


0.3510 


0.6150 


0.4224 


Oil 


0.6864 


0.4689 


0.4718 


0.6983 


0.4710 


0.4704 



Table 2 Atomic positions for the high and low temperature structure of VeOn as used in the band structure calcula- 
tions. These data have been determined by Canfield [60]; the coordinates refer to the primitive translations proposed 
by Horiuchi et al. [62]. 



The low temperature structure at 20 K results in (A = 8.6981 as) 

/ -1.0000 \ / 1.0000 \ / -0.0074 \ 

a H = A I 0.0000 , b H = A I 0.9224 , c H = A I 5.3417 . (8) 
V 0.6521 J \ 0.6521 / \ 3.7636 / 

The Bohr radius is denoted by as and the primitive translations refer to Cartesian coordinates. The corre- 
sponding positional parameters of the vanadium and oxygen atoms are denoted in table|2] 

The presented LDA calculations are based on the (scalar relativistic) augmented spherical wave (ASW) 
method [66, 67]. We adopt the parametrization of the exchange correlation potential introduced by Vosko, 
Wilk, and Nusair [68]. The ASW scheme uses the atomic sphere approximation and models the full crystal 
potential by means of spherical symmetric potential wells. Here it is required that the atomic spheres fill the 
space of the unit cell. For open crystal structures problems can arise because space filling only due to atom 
centered spheres may lead to large overlap. As a consequence, so-called empty spheres, i.e. pseudo atoms 
without nuclei, have to be introduced. They are used to correctly model the shape of the crystal potential in 
large voids. The collection of both physical and empty spheres leads to an artificial close-packed structure. 
However, since the potential of the whole set of spheres ought to represent the crystal potential as exactly 
as possible it is a challenge to find optimal empty sphere positions and optimal radii for the real and empty 
spheres. Here the sphere geometry optimization algorithm described in [69] is most efficient. 

By adding 35 empty spheres from 18 crystallographically inequivalent classes to the triclinic unit cell of 
high temperature VgOn it is possible to keep the linear overlap of real spheres below 18%. Simultaneously, 
the overlap of any pair of real and empty spheres is smaller than 23%. In the case of the low temperature 
structure 40 empty spheres from 21 inequivalent classes allow for reducing the overlaps to less than 18% 
and 23%, respectively. Summing up, the unit cell entering the LDA calculation comprises 69 spheres for 
the high and 74 spheres in the case of the low temperature structure. The radii of the vanadium and oxygen 
spheres are summarized in table |3j which additionally denotes the valence charges arising from the LDA 
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High temperature structure 


Low temperature structure 


Atom 


Radius 


Charge 


Radius 


Charge 


VI 


2.1866 


2.5090 


2.2303 


2.6179 


V2 


2.1587 


2.7155 


2.2278 


2.8276 


V3 


2.2833 


2.5240 


2.2158 


2.5534 


V4 


2.2755 


2.7564 


2.2248 


2.7939 


V5 


2.3440 


2.7775 


2.2248 


2.5213 


V6 


2.3437 


2.7354 


2.1856 


2.6201 


Ol 


1.9684 


4.1741 


1.9423 


4.2072 


02 


1.7883 


4.0164 


1.8240 


3.9466 


03 


1.9680 


4.3682 


1.9699 


4.4187 


04 


2.0093 


4.3624 


1.9420 


4.2381 


05 


1.7654 


3.8853 


1.8219 


3.8832 


06 


1.8752 


4.0164 


1.8201 


3.9683 


07 


1.8673 


4.0673 


1.7943 


3.9711 


08 


1.8908 


4.0278 


1.8394 


4.0294 


09 


1.8585 


4.0833 


1.8568 


4.0498 


O10 


1.8828 


4.1106 


1.8121 


3.9479 


Oil 


1.9170 


4.0959 


1.8166 


3.8850 



Table 3 Radii of the vanadium/oxygen spheres (in as) as well as calculated LDA valence charges (V 3d or O 2p) for 
both the high and the low temperature phase of VgOn. 



band structure calculation. For both VgOn structures the basis sets taken into account in the secular matrix 
comprise V 4s, Ap, 3d, (4/) and O 2s, 2p, (3d) orbitals. States given in parentheses enter as tails of other 
states; for details on the ASW method see [67,70]. To complete the basis sets we add empty sphere states, 
which are determined with respect to the spacial extensions of the spheres and the total charge occupying 
them. The used configurations reach from Is, (2p) to Is, 2p, 3d, 4/, (5g). During the course of the LDA 
calculation the Brillouin zone is sampled with an increasing number of k-points in the irreducible wedge. 
In this manner one ensures convergence of the results with respect to the fineness of the k-space grid. For 
both the high and low temperature calculation the number of k-points was increased from 108 to 256, 864, 
and 2048. Self-consistency of the charge density was assumed for deviations of the atomic charges and the 
total energy of subsequent iterations less than 10 -8 electrons and 10~ 8 Ryd, respectively. 

As V4O7 reveals short metal chains (length n = 4) it is suitable for investigating the electronic features 
of the sesquioxide-like chain end sites and their influence on the MIT. Using our results for VgOn we can 
understand the MIT of V4O7. Afterwards we establish important implications for the frequently discussed 
MITs of V 2 3 . As for V 6 On, the space group of V4O7 is PI {C}) and the MIT at 250 K is accompanied 
by a structural distortion. Refinements tracing back to Hodeau and Marezio [59] placed the vanadium and 
oxygen sites at the Wyckoff positions (2i): ±(x, y, z). Above the transition at 298 K the authors observed 
the triclinic parameters a A = 5.509 A, b A = 7.008 A, c A = 12.256 A, a A = 95.10°, (3 A = 95.17°, and 
7a = 109.25°. The corresponding values at low temperatures (120 K) are a A — 5.503 A, b A = 6.997 A, 
c A = 12.256 A, a A = 94.86°, (i A = 95.17°, and j A = 109.39°. This refers to the unit cell of Andersson 
and Jahnberg, see equation Q with n = 4. The primitive translations of Horiuchi et al. are given by 



SLH \ ( -1 1 \ / Blr \ ( 1 \ / 8LA \ 

b H = 1 1 1 b fl = 1 b A 

oh J V I I J \ or J V 2 § 4 / \ c J 
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High temperature structure 


Low temperature structure 


Atom 


X 


y 


z 


X 


V 

a 


z 


VI 


0.9806 


0.9383 


0.1008 


0.9773 


0.9299 


0.1043 


V2 


0.9946 


0.4509 


0.0964 


0.9781 


0.4382 


0.1008 


V3 


0.9662 


0.9613 


0.3742 


0.9508 


0.9512 


0.3757 


V4 


0.9914 


0.4905 


0.3658 


0.0146 


0.5021 


0.3608 


Ol 


0.6927 


0.4137 


0.0514 


0.6932 


0.4151 


0.0512 


02 


0.3237 


0.5551 


0.1036 


0.3191 


0.5383 


0.1112 


03 


0.6825 


0.4643 


0.1716 


0.6808 


0.4658 


0.1718 


04 


0.3763 


0.6232 


0.2254 


0.3804 


0.6214 


0.2254 


05 


0.6893 


0.4112 


0.3328 


0.7018 


0.4240 


0.3278 


06 


0.3044 


0.5683 


0.3910 


0.2793 


0.5257 


0.4068 


07 


0.6605 


0.4326 


0.4700 


0.6616 


0.4218 


0.4728 



Table 4 Atomic positions for the high and low temperature structure of V4O7 as used in the band structure calcula- 
tions. These data trace back to Hodeau and Marezio [59]; the coordinates refer to the primitive translations proposed 
by Horiuchi et al. [62]. 



and we are able to transform coordinates from the A-system into the /J-system 

(10) 




Again, the translation vector has been added to account for different choices of the origin in the coordinate 
systems. Using equation d 1 01 the positional parameters of Hodeau and Marezio for high and low temper- 
ature V4O7 result in the values summarized in table |4] which enter the LDA calculation. Alternative data 
for the crystal structure were given by Horiuchi et al. [57] and Marezio et al. [58]. Altogether, one obtains 
four crystallographically inequivalent vanadium and seven oxygen sites. Due to the inversion symmetry of 
the lattice the V4O7 unit cell contains 22 atoms. 

Using the lattice constants of Hodeau and Marezio we set up the triclinic unit cell of V4O7 in Cartesian 
coordinates. Because of the similar definitions of the A and L-system we can proceed analogous with our 
considerations for VgOn, applying the parameters a a, i>A, ca, &a, Pa and ^a- Via equation (|9jl we find a 
unit cell in the iJ-system. For the high temperature structure at 298 K this yields {A = 8.7702 as) 

-1.0000 \ / 1.0000 \ / 0.0071 

& H = A\ 0.0000 , b H = A\ 0.9333 , c H = A [ 3.4280 | . (11) 
0.6396 / \ 0.6396 / \ 2.3282 

Moreover, the low temperature configuration at 120 K results in (A = 8.7691 as) 

I -1.0000 \ / 1.0000 \ / 0.0054 

sljj = A I 0.0000 , b H = A\ 0.9313 , c H = A [ 3.4199 | . (12) 
\ 0.6374 / \ 0.6374 / \ 2.3141 

These primitive translations are used with the positional parameters denoted in table|4] 

Similar to the procedure in the case of VgOn we insert additional augmentation spheres into the V4O7 
structure to account for its openness and to properly model the crystal potential. For both the high and the 
low temperature phase it suffices to apply 13 crystallographically inequivalent classes and place altogether 
24 empty spheres in the triclinic unit cell. This keeps the linear overlap of physical spheres below 19% and 
the overlap of any pair of physical and empty spheres below 24%. Unit cells entering the LDA calculation 
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High temperature structure 


Low temperature structure 


Atom 


Radius 


Charge 


Radius 


Charge 


VI 


2.2002 


2.7937 


2.1388 


2.7457 


V2 


2.2349 


2.6572 


2.3468 


2.7529 


V3 


2.3224 


2.7329 


2.1999 


2.7430 


V4 


2.3433 


2.7245 


2.3436 


2.6778 


Ol 


1.9931 


4.3397 


2.0096 


4.3259 


02 


1.8278 


4.0266 


1.9267 


4.1770 


03 


1.9883 


4.4309 


1.9759 


4.4115 


04 


2.0465 


4.3530 


2.0226 


4.3440 


05 


1.7994 


3.9182 


1.7492 


3.8328 


06 


1.8993 


4.0359 


1.8357 


3.9669 


07 


1.9164 


4.1097 


1.9045 


4.0616 



Table 5 Radii of the vanadium/oxygen spheres (in as) as well as calculated LDA valence charges (V 3d or O 2p) for 
both the high and the low temperature phase of V4O7. 



comprise 46 spheres; table |5]denotes the radii of the physical spheres and the calculated valence charges. 
All technical details of the calculation are the same as for VgOn. 

By comparing V7O13 to VgOis the influence of the chain center atoms on the electronic properties of 
the Magneli phases will be investigated quantitatively in section0 Such a comparison is of special interest 
since the former compound is the only member of the Magneli series, which does not undergo an MIT. In 
contrast, V 8 Oi5 exhibits an MIT at » 70 K accompanied by structural distortions. V7O13 and both VgOis 
phases crystallize in the triclinic space group PI (C]) with the vanadium and oxygen atoms at the Wyckoff 
positions (2i): ±(x, y, z). Canfield [60] reported lattice parameters related to the unit cell of Le Page and 
Strobel, see equation (0. For V 7 Oi 3 wehavea L = 5.439 A, b L = 7.013A, c L = 37.161 A, a L = 67.04°, 
Pl = 57.46°, and 7^ = 108.92°. For the room temperature structure of V 8 Oi5 these constants amount to 
a L = 5.431 A, b L = 7.017 A, c L = 42.896 A, a L = 66.84°, (3 L = 57.55°, and j L = 108.94°, whereas 
the low temperature configuration is defined by the values — 10.892 A, — 6.980 A, cl = 85.907 A, 
qj, = 66.83°, (3l = 57.38°, and 7^ = 108.67°. The latter unit cell is exceptional because it accounts for 
a superstructure evolving in V 8 Oi5 at l° w temperatures, which complicates the considerations. However, 
for the high temperature structures we proceed analogous with the previous analysis of VgOn. A detailed 
discussion of the structural input for the subsequently presented electronic structure calculations has been 
given in [7 1 ] . We use the positional parameters measured by Canfield. Alternative cry stallographic studies 
of V7O13 at room temperature were performed by Horiuchi et al. [62]. 

5 Dimerization and localization in V 6 On 

In this section we analyze the MIT of VgOn [72]. Changes of the electronic structure at the transition are 
discussed in relation to the structural transformations occuring simultaneously. The study will benefit from 
our unified representation of the crystal structures of the Magneli phases as well as of VO2 and V2O3. We 
will succeed in grouping the electronic bands of V 6 On into states behaving similarly to either the dioxide 
or the sesquioxide. Therefore it is possible to analyze the phase transitions of the different compounds on 
the basis of a common point of view, which helps us to gain insight into the delicate interplay of electron- 
lattice coupling and electronic correlations and their influence on the respective MITs. 

We start our considerations with a survey of the partial densities of states shown in figure^] which arise 
from the LDA calculations for high and low temperature VgOn. The low temperature antiferromagnetism 
is not considered for the reasons discussed earlier. In the high temperature case we observe three groups of 
bands in the energy intervals from —8.2 eV to —2.9 eV, from — 0.7eV to 1.8 eV, and from 2.0 eV to 4.8 eV. 
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Fig. 10 Partial V 3d and O 2p densities of states (DOS) per vanadium atom resulting from the high (298 K) and low 
(20 K) temperature crystal structure of VeOn are compared. 



On entering the low temperature phase the lowest group broadens considerably, reaching from — 8.4 eV to 
— 2.7eV, whereas the other groups are almost not altered. Due to substantial structural similarities to the 
dioxide/sesquioxide it is not surprising that the electronic structure of VeOn resembles the expectations of 
the molecular orbital picture as discussed for the former compounds. Characteristic a and 7r-type overlap 
of V 3d and O 2p states gives rise to bonding-antibonding split molecular orbitals in the energetical order 
f7-7r-7r*-(T*. Bonding states are identified with the lowest, tt* states with the middle, and a* states with the 
highest group of bands. Since the metal atoms are octahedrally coordinated, the a* and 7r* bands display 
e° and t^g symmetry, respectively. Because of 12 V and 22 O atoms per unit cell there are 12 x 2 — 24 
vanadium 3d & a g , 12 x 3 = 36 vanadium 3d t2 g , and 22 x 3 = 66 oxygen 2p bands. In agreement with the 
molecular orbital picture and the findings for VO2 and V2O3 the energetically lowest structure in figurelTUl 
predominantly traces back to O 2p states, whereas the other groups mainly originate from V 3d states. Due 
to a different electron count, the filling of the t^g levels is between that of VO2 and V2O3. Contributions of 
vanadium and oxygen in regions dominated by the respective other states are due to p-d hybridization. In 
regions corresponding to er-type V-O overlap, i.e. in the e g range and in the lower half of the O 2p range, 
they are stronger. Larger t2 g admixtures in the e g interval of VgOn, compared to VO2, point to increased 
deviations from an ideal octahedral coordination. The distortions agree more with those of V2O3. 

Differences between the partial densities of states calculated for the high and low temperature structure 
are small. No energy gap is observed for the low temperature structure. Instead, the structural modifications 
accompanying the MIT leave the DOS at the Fermi energy almost unaffected. This is not surprising as for 
neither the insulating phase of VO2 nor the AFI phase of V2O3 did the calculations succeed in reproducing 
an energy gap, which fact was attributed to shortcomings of the LDA. However, as has been demonstrated 
for the dioxide, this does not prevent us from understandig the mechanism of the MIT [24, 36]. Since we 
study the relations between electronic states and local atomic environments as well as the modification of 
this relationship at the transition, the LDA limitations do not affect the following considerations. 

The structural modifications of VgOn at 170 K are mainly concerned with a strong V4-V6 dimerization 
in the 2-4-6 vanadium chain, see figure|5] While the V4-V6 bond length decreases from 2.85 A to 2.62 A, 
the V2-V4 and V6-V6 distances grow from 2.95 A to 3.06 A and from 2.82 A to 3.30 A, respectively. As a 
consequence, two isolated V4-V6 vanadium pairs arise and the chain end atoms (V2) separate from the rest 
of the vanadium chain. In the 1-3-5 chain dimerization effects yield a shortened V5-V5 bond length both 
above and below the MIT. Obviously, in the insulating phase the V5-V5 pair separates even more from the 
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Fig. 11 Site-projected partial V 3d t2 g densities of states (DOS) per metal atom for the high and the low temperature 
crystal structure of VgOii: sites VI and V2. The crystal structure is depicted in figure|9|and the orbitals refer to the 
local rotated reference frame. 



evolving VI -V3 pairs. All these modifications of the crystal structure resemble the pairing effects known 
from the Peierls distortion in monoclinic VO2. Due to finite chains in the case of VgOn it is not surprising 
to find more complicated distortion patterns; the end atoms affect the V-V pairing in the whole chain. 

The monoclinic phase of VO2 is characterized not only by metal-metal pairing but also by zigzag-type 
in-plane displacements of the vanadium atoms along the local z-axis, thus along the diagonals of the rutile 
basal planes. While metal sites in the rutile structure coincide with the centers of the surrounding oxygen 
octahedra, they are shifted away from these positions by 0.20 A in the distorted configuration, thus giving 
rise to stronger V-O overlap. Equivalent shifts away from the centers of the octahedra are observed for the 
sites V3, V4, V5, and V6 in VgOn. The distance between vanadium atoms and octahedral centers grows 
from 0.16 A to 0.29 A (V3), from 0.17 A to 0.23 A (V4), from 0.07A to 0.23 A (V5), and from 0.04 A to 
0.35 A (V6), respectively. Actually, the major part of the shifts is oriented perpendicular to the c pm t-axis. 
Although the metal sites in the high temperature phase do not coincide with octhedral centers, the relative 
changes at the MIT are similar to VO2. Interestingly, the high temperature lateral displacements are larger 
for the intermediate sites V3 and V4 than for the center atoms V5 and V6. This is due to a slight rotation 
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Fig. 12 Same representation as used in figure fTTI but for the vanadium sites V3/V4. 



of the vanadium chains away from the c pru t-axis, which is present in each Magneli phase. Because of the 
rotation the chain end sites VI and V2 shift away from their nearest neighbours along a prut - resembling 
the vanadium anti-dimerization of V2O3 parallel to the Chex-axis. Here the displacements of the metal sites 
with respect to octahedral centers amount to 0.18 A (PM phase) and 0.21 A (AFI phase). For VgOn these 
values are 0.32 A in the metallic and 0.35 A in the insulating phase. Hence, even concerning the details of 
the local distortion, chain centers and ends behave as their parent structures VO2 and V2O3. 

Next we decompose the V 3d t2 g group of states into its symmetry components. FiguresfTT1toll3lgive 
site-projected partial d x 2_ y 2, d yz , and d xz densities of states for the six inequivalent metal sites V1,...,V6. 
For each site the presentation of the data refers to the local rotated reference frame, which was defined in 
the discussion of the dioxide. The rutile reference frame is useful since the local octahedral coordinations 
of the metal sites in VgOn resemble the rutile arrangement, see the preceding sections. 

We turn to the <2g states of the vanadium atoms V4, V5, and V6, which are involved in strong dioxide- 
like displacements. Thus the site-projected densities of states are similar to one another and resemble the 
DOS of rutile and monoclinic VO2. Note the two peak structure weakly indicated in the high temperature 
V4 d x 2_ y 2 DOS due to cr-type V-V bonding along the metal chains. In the low temperature phase the V4- 
V6 dimerization causes an increased splitting of the DOS into bonding and antibonding branches located 
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Fig. 13 Same representation as used in figure fTTI but for the vanadium sites V5/V6. 



at energies of — 0.5eV and 1.4eV. The d yz and d xz densities of states undergo energetical upshifts due to 
the antiferroelectric displacements of the V4 atoms perpendicular to c prut , which raise the overlap of V 3d 
and O 2p states and thus the 7T-7T* splitting. Similar to VO2, the energetical separation of the d x 2_ y 2 = d\\ 
band and the remaining e£ states is increased but not complete in the low temperature results. 

The situation is similar for the V6 atoms, which participate in the V4-V6 pairs. Therefore the bonding 
and antibonding branches of the d x i_ y i DOS appear at the same energies as for the V4 site. However, in 
the low temperature configuration the splitting is more pronounced due to additional V6-V6 bonding. As 
in the case of V4 the energetical upshift of the d yz ld xz DOS is easily observed, especially directly below 
the Fermi level. While the d xz DOS of NAN 6 gives rise to a broad single peak, we find a distinct splitting 
of the d yz DOS in contributions centered at about — 0.2 eV and 1.0 eV for the high temperature structure. 
The splitting is understood by taking into account metal-metal bonding parallel to 6 pru t. Across octahedral 
faces typical sequences V2-V6-V4 appear in this direction, as obvious from figure[6] Due to a peak in the 
low temperature V6 d yz DOS at about 1.4 eV, which is not present in the respective V4 DOS but reappears 
in the V2 DOS, we can assume increased V2-V6 bonding along 6 pru t at low temperatures. The d x 2_ y 2 and 
d yz densities of states of atom V5 show strong bonding-antibonding splitting, which hardly changes at the 
transition. The shape of the d x 2_ y 2 DOS is understood in terms of constant V5-V5 bonding. In contrast, 
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splitting of the d yz DOS traces back to V-V overlap along & prut . The energetical upshift of these states at 
the MIT is pronounced and has the same origin as discussed for V4 and V6. In total, the partial V4, V5, 
and V6 densities of states display characteristic features known from the dimerization and the zigzag-type 
displacements in monoclinic VO2 and are understood in terms of the latter compound. 

In contrast to the chain center atoms V4, V5, and V6 the end atoms VI and V2 are characterized by a 
sesquioxide-like local environment. Thus they are involved in V-V bonding across shared octahedral faces 
connecting the 1-3-5 to the 2-4-6 layers, but they are not subject to dimerization along c prut . Hence, except 
for small peaks and shoulders around — 0.5eV and from 1.0 eV to 1.3 eV (reminiscent of the chain center 
d x 2_ y 2 densities of states) the V1/V2 d x 2_ y 2 DOS consists of a single broad peak extending from — 0.2 eV 
to 0.7 eV. The subpeaks observed in this energy region miss counterparts in the DOS of any neighbouring 
vanadium atom. Since such peaks are indicative of metal-metal bonding we conclude, in consistence with 
the reported V-V distances, that chain end atoms are well separated from the remainder of the chains. Thus 
the d x 2_ y 2 states can be regarded as localized. While chain center d x 2_ y2 orbitals predominantly mediate 
V-V overlap along c prut , this interaction is suppressed for VI and V2 for geometrical reasons. Hence the 
chain end d x 2_ y 2 states behave as the d xz states, which generally give rise to reduced (7r-type) metal-metal 
overlap along the chains. In both phases the center of weight of the VI d x 2_ y 2 DOS is noticeably higher 
than that of the V2 d x 2_ y 2 DOS. Because of nearest neighbours along a pru t the d yz orbitals are important 
for interpreting the chain end atoms since these states cause V-V bonding parallel to both a prut and 6 pru t- 
The VI and V2 d yz states display bonding-antibonding splitting both above and below the transition. For 
the insulating modification this separation is more pronounced with peaks directly above the Fermi energy 
and around 1.2 eV. However, so far we are not able to distinguish between the effects of the different metal- 
metal interactions. Whether the shape of the d yz DOS indicates bonding along a prut or 6 prut or both will 
be dealt with later. The question is of interest as the a prut -direction of V 6 On corresponds to the q lox -axis 
of the corundum structure. Hence the splitting of the d yz states is equivalent to the a\ g splitting in V2O3 
and we can gain new insight into the MIT of the sesquioxide. The partial VI and V2 t 2g densities of states 
resemble each other due to similar coordinations. They and are interpreted in accordance with V2O3. 

Atom V3 is exceptional as its t 2g DOS can neither be interpreted by means of dioxide-like dimerization 
in the vanadium chains nor by sesquioxide-like metal-metal overlap along a prut . The second fact is simply 
due to the absence of adjacent vanadium sites, which allow for overlap between the 1-3-5 and 2-4-6 layers. 
Nevertheless, all components of the t 2g DOS are similar to those of sites VI and V2. The d x 2_ y2 and d xz 
densities of states consist of broad single peaks. Reminiscent of bonding and antibonding branches in the 
corresponding V5 DOS we find satellite peaks at — 0.4eV and a shoulder at 1.0 eV in the high temperature 
phase. The compact shape of the d xz DOS again is not surprising due to the small V-V overlap this orbital is 
involved in, not only for V3 but also for the other vanadium atoms. In contrast, the d yz partial DOS reveals 
a distinct double peak structure in both phases, indicative of bonding-antibonding splitting. At first glance 
this finding is difficult to understand as geometrically no vanadium partner along a prut exists, which could 
participate in V-V bonding. However, this puzzling situation is resolved by reconsidering the possibility of 
V-V interaction perpendicular to the chains and simultaneously perpendicular to the ai g -direction, namely 
along 6 prut . As mentioned before, there are finite vanadium sequences consisting of three atoms (V1-V5- 
V3), where the distance between neighbouring atoms corresponds to the longer rutile lattice constant. To 
be more precise, these bond lengths amount to 4.58 A/4. 54A (V1-V5) and 4.45 A/4.40 A (V5-V3) for the 
high/low temperature configuration. The associated values in the V2-V6-V4 sequences are 4.54 A/4.57 A 
(V2-V6) and 4.49 A/4.38 A (V6-V4). Since each vanadium atom in VeOn participates in er-type overlap 
along fe pru t, the bonding-antibonding splitting of the d yz states is present not only for site V3. 

This has important implications for the interpretation of the V1/V2 d yz DOS because interactions par- 
allel to a pru t and 6 pru t add to the bonding-antibonding splitting of these states. The respective antibonding 
peaks of VI and V2 are merged at high temperatures but split up at the MIT. Due to the V-V overlap along 
frprut site V3 resembles the electronic properties of the chain end sites, although its coordination does not 
allow for ai g -like bonding. Therefore V-V bonding across the vanadium layers is not the only source for 
a splitting of the d yz DOS. Despite differing bond lengths hardly any difference between V-V overlap via 
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octahedral faces within and perpendicular to the layers is found. Supporting results of Elfimov et al. [52], 
V-V couplings other than the a prut -like are important for the shape of the a\ g bands in V2O3. 

While V-O overlap places the V 3d ti g states near the Fermi level, the detailed electronic features of the 
orbitals, and thus the MIT of VgOn, are fundamentally influenced by the local metal-metal coordination. 
The small variations of the oxygen sublattice at the MIT barely affect the shape of any t 2g DOS, whereas 
modifications of the V-V bond lengths are important. Depending on the coordination, the local electronic 
properties of the vanadium atoms are typical for VO2 or V2O3 and it is reasonable to regard the partial den- 
sities of states as local quantities. At the chain centers (V4, V5, V6) dimerization and antiferroelectric-like 
displacements via strong electron-lattice interaction cause a splitting of the d x 2_ y 2 DOS and an energetical 
upshift of the e g states at the MIT. For atom V5 the former effect is already present at high temperatures. 
Summarizing, chain center atoms behave analogous with VO2. In contrast, at the chain ends (VI, V2, V3) 
the d x 2_ y 2 orbitals are characterized by strongly reduced metal-metal overlap and show a localized nature. 
The shape of the d x 2_ y 2 DOS is closely related to the d xz DOS. Due to the strong localization, electronic 
correlations are important for the chain end d x 2_ y 2 states. Being subject only to in-layer V-V overlap, the 
partial V3 d yz DOS shows bonding-antibonding splitting similar to the atoms VI and V2, which have an 
additional nearest neighbour across the layers. Therefore the in-plane V-V interaction is at least as impor- 
tant as the perpendicular overlap. In conclusion, the sesquioxide-like regions of the VgOn crystal appear 
to be susceptible to electronic correlations for geometrical reasons. Structural variations at the MIT leave 
all three components of the t 2g partial DOS almost unchanged, thus reflecting close relations to V2O3. We 
interpret the phase transition of VgOn as resulting from a combination of electron-lattice interaction and 
electronic correlations. Due to the different atomic arrangements in VO2 and V203-like regions neither an 
embedded Peierls instability nor correlations can account for the MIT of VgOn on their own. 

Going back to the calculated valence charges at the vanadium and oxygen sites summarized in table[3]we 
investigate possible charge ordering. Different atomic spheres radii, necessary to fulfill the requirements of 
the atomic sphere approximation, prohibit a quantitative comparison of the valence charges but allow still 
for qualitative results. For the high temparature modification there are about 0.2 electrons less at sites VI 
and V3 than at sites V2 and V4, whereas the charges at the central sites V5 and V6 coincide. The atomic 
spheres comprise almost the same space for the pairs V1/V2, V3/V4, and V5/V6, respectively. In general, 
the calculated 2.5/2.7 V 3d valence electrons contradict the simple ionic picture predicting an average of 
1.33 electrons per vanadium atom. Instead of V 3 67+ ions we find V 2 3+ and V 2 - 5+ configurations, which 
is not surprising due to covalent V-O bonding. Oxygen spheres comprise about four 2p electrons. Further 
electrons enter the empty spheres and hence the open space between vanadium and oxygen. 

The accuracy of comparing calculated valences is limited as the assignment of charge to specific sites is 
arbitrary to some degree. We assign charge within an atomic sphere to the respective atomic site. Electrons 
entering interstitial empty spheres cannot be assigned to any specific atom. Importantly, differences in the 
occupation of equivalent states can be identified with a much higher accuracy. Thus the 3d charge of both 
VI and V3 is significantly reduced. In addition, we obtain the same result for the low temperature phase, 
where all vanadium spheres show similar radii. The valence charge at sites VI, V3, and V5 is roughly 0.2 
electrons smaller than at sites V2, V4, and V6, resulting in a charge transfer of 0.1 electrons per site from 
the 1-3-5 to the 2-4-6 chain above and below the MIT. The different electron count in the two metal chains 
is related to their crystallographic inequivalence. While donor atoms (to some extend) approach the VO2 
charge configuration, the acceptor chain resembles more the V2O3 filling. The electronic configurations of 
VO2 and V2O3 seem preferable to a mixture of both. However, chain end sites do not generally gain excess 
charge at the expense of the dioxide-like chain centers. Hence it is not the electron count of a specific atom 
but the local coordination which determines the behaviour at the phase transition. The calculated valence 
charges contradict V 3+ /V 4+ states, but confirm two kinds of metal valences; they point to smaller charge 
differences. Deviations from the ionic picture are consistent with susceptibility data by Gossard et al. [13], 
whereas for the insulating phase we do not find the proposed charge ordering. 
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The above findings have important implications for understanding V2O3 as the chain end atoms of the 
Magneli phases resemble the local metal coordination of this material. According to the previous discus- 
sion the electronic properties of VgOn are influenced by essentially three effects. First, the center atoms 
of the vanadium chains along c pru t reflect characteristics of the embedded Peierls instability responsible 
for the MIT in V02- Second, the localization of the chain end d x 2_ y 2 orbitals leaves them susceptible to 
electronic correlations. Third, the chain end d yz orbitals to roughly equal parts are subject to both metal- 
metal bonding within (along 6 pru t) and perpendicular to (along a pmt ) the vanadium layers. As vanadium 
chains in the case of V2O3 (n — 2) degenerate to pairs, a dioxide-like Peierls instability cannot contribute 
as a matter of principle. However, the localized nature of the d x i_ y i orbitals and the importance of strong 
V-V coupling along & pm t should be transferred to the sesquioxide. Due to the high crystal symmetry both 
effects could not be analyzed in previous studies of V2O3. As a consequence of the broken symmetry in 
the Magneli phases, combined with a closely related coordination of appropriate vanadium atoms, they are 
accessible to the study of VgOn. The splitting of the ai 9 -like states due to metal-metal overlap along Chcx 
therefore seems to be less important than commonly assumed for V2O3. 

6 Band narrowing in V4O7 

The compound V4O7 is suitable for studying chain end effects since the metal chains here comprise only 
four atoms [73]. Analogous with the preceding section we discuss the electronic structure and its changes 
at the MIT in relation to the simultaneous structural transformations. Applying the knowledge established 
in the study of VgOn we can group the electronic bands into states behaving either VO2 or V203-like. 

Partial V 3d and O 2p densities of states arising from the high and low temperature structure are shown 
in figure^] As for VO2, V2O3, and VgOn, we find three separated groups of electronic bands. Recalling 
the molecular orbital picture, the bonding-antibonding splitting due to V-O overlap, and the discussions in 
the previous sections, we identify the bonding states with those bands extending from — 8.3 eV to — 3.3 eV 
at high temperatures; they shift to the region from — 8.2eV to — 3.0eV below the MIT. Further groups of 
bands correspond to ir* and a* states. Hence they are dominated by vanadium states and separate into t2 g 
and gg contributions due to the octahedral coordination of the metal sites. While the t2 g group is observed 
between — 0.7eV and 1.6eV in the metallic phase, its upper edge shifts to 1.7eV at low temperatures. The 
e y group starts at 2.0 eV in both phases, but it extends to 4.7 eV instead of 4.4 eV in the insulating phase. 
Since the unit cell of V4O7 contains 8 vanadium and 14 oxygen atoms we have 8x2 = 16 vanadium 3d 
e^, 8 x 3 = 24 vanadium 3d ti g , and 14 x 3 = 42 oxygen 2p bands. Hybridization between V 3d and 
O 2p orbitals results in finite contributions of vanadium in the bonding and of oxygen in the antibonding 
region. As expected, such hybridization effects are stronger for er-type V-O overlap. Moreover, there is a 
finite <2g DOS in the e a g dominated region. As the magnitude of this admixture equals the result for VgOii, 
deviations from an ideal octahedral coordination of the vanadium sites are similar in both compounds and 
likewise similar to the sesquioxide. All these results for V4O7 fit very well into the picture developed in 
the preceding sections. Due to the formal V 3d charge of 1.5 electrons per metal site, V4O7 systematically 
is between VgOn and V2O3. Therefore the V 3d t^g orbitals are filled to a somewhat larger extent than 
for VgOn, which is confirmed by fi gures fTol and 1 1 41 Although the calculation misses an energy gap for 
the experimental insulating phase, the LDA shortcomings again do not affect an investigation of relations 
between structural and electronic properties or the analysis of implications for the MIT. 

While the previous analysis of VgOn gave strong indications for a predominant influence of electron- 
lattice interaction in dioxide-like and of electronic correlations in sesquioxide-like regions of the crystal, 
we expect that the investigation of the sesquioxide-related member V4O7 will allow for conclusions about 
the MIT of V2O3 . Separated by oxygen atoms, inequivalent vanadium layers alternate along the a pru t-axis 
of V4O7, comprising the atoms V1/V3 or V2/V4. Due to this alternation and due to relative shifts of the 
metal 4-chains along c prut the end sites VI and V2 are located on top of each other, see figure|8] We again 
refer the V 3d orbitals to local coordinate systems with the z-axis oriented parallel to the apical axis of the 
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Fig. 14 Partial V 3d and O 2p densities of states (DOS) per vanadium atom resulting from the high (298 K) and low 
(120 K) temperature structure of the compound V4O7. 



local octahedron and the x-axis parallel to c pru t, respectively. Using this reference frame, a decomposition 
of the V 3d ti g group of states into the symmetry components d x 2_ y 2, d yz , and d xz is depicted in figures 
E]and^] The site-projected partial densities of states are given for the vanadium sites VI, V2, V3, and 
V4. In the following, we concentrate on the results for the 1-3 chain, which is motivated by the fact that 
the modifications of the local surroundings of atom VI at the MIT are similar to those known from V2O3. 
According to the data denoted in figure[8]the V1-V2 and V1-V3 distances increase at the MIT, as do the 
distances a and b in the sesquioxide at the PM-AFI transition, see figure|4] Since in the latter compound all 
vanadium atoms are crystallographically equivalent a complete symmetry analysis of the electronic states 
cannot be performed. However, this is possible for V4O7. As demonstrated in the case of VgOn, the local 
atomic arrangements are related to the electronic structures of the respective sites. Therefore the VI site of 
V4O7 allows us to study the influence of the structural changes at the MIT of V2O3 in more detail than in 
a direct investigation. 

Changes in the V4O7 crystal structure at 250 K are mainly concerned with strong dimerization evolving 
in both the 1-3 and 2-4 metal chain, see figure[8] While the V3-V3 and V2-V4 bond lengths shorten from 
2.79 A and 2.93 A to 2.67 A and 2.83 A, the remaining distances V1-V3 and V4-V4 elongate from 2.97A 
and 2.93 A to 3.03 A, respectively. Hence we find an isolated V3-V3 pair in the first and two V2-V4 pairs 
in the second chain. In the 1-3 chain the end atoms decouple to a large degree from the central V3-V3 pair. 
The modifications resemble pairing effects associated with the embedded Peierls instability of the dioxide. 
Nevertheless, the end atoms of the 1-3 chain likewise show structural changes known from the symmetry 
breaking PM-AFI transition of V2O3. The distance between VI and the nearest vanadium site along c pru t 
increases by 0.06 A, which is half the shift found for the sesquioxide. Note that the dominating structural 
changes in V2O3 do not appear along a pru t but along c pru t. In contrast to the 1-3 chain, the end atoms of 
the 2-4 chain increase the coupling to adjacent in-chain sites below the MIT. Thus they show the opposite 
trend than found in V2O3 and are not suitable for investigating the latter material. 

Beyond metal-metal pairing also zigzag-type in-plane displacements of the metal sites are present in the 
V4O7 structure. The distances between vanadium atoms and octahedral centers amount to 0.31 A/0.35 A 
(VI), 0.32A/0.27A (V2), 0.15 A/0.17A (V3), and 0.08 A/0.09 A (V4) in the high/low temperature modi- 
fication. Apparently, the average bond length does not increase at the transition, which contrasts findings 
for VO2 and VgOn, hence confirming the sesquioxide-like character of V4O7. Because of a slight rotation 
of the vanadium chains away from the c prut -axis, the lateral displacement is much larger for the chain end 



adp header will be provided by the publisher 



27 



DOS ( 1/eV) 
3 



298 K 



V 4 7 



120 K 



VI site 


V 3d K ! ; 

V3d„ ■ 

V3d„ " 

















DOS (1/eV) 
3 



1 
(E - E F ) (eV) 

298 K 



V 4 7 



1 2 

(E - E F ) (eV) 



120 K 



V2 site 


V 3d x 2. y 2 

V3d y2 ■ 

V3d„ ■ 















I 

(E - E F ) (eV) 



1 

(E - E F ) (eV) 



Fig. 15 Site-projected partial V 3d t2 g densities of states (DOS) per metal atom for the high and the low temperature 
crystal structure of V4O7: sites VI and V2. The crystal structure is depicted in figure[8|and the orbitals refer to the 
local rotated reference frame. 



sites VI and V2. The rotation enables the latter sites to increase their distance along a pru t, which reflects 
the anti-dimerization of V2O3 parallel to Ch ex - In general, the main part of the zigzag-type displacements is 
oriented perpendicular to c prut . Atom V2 shifts almost along the apical axis of its octahedron, i.e. along the 
local z-axis, whereas atom VI moves perpendicular to the local z-axis. However, because of the rotation 
of the local coordinate system by 90° around the c prut -axis, the displacements occur in every chain along 
the same direction. On closer inspection we recognize the motions to be antiparallel, which increases the 
in-pair VI -V2 bond length. Altogether, the zigzag-type antiferroelectric mode of the 2-4 chain resembles 
the behaviour of the dioxide, whereas the 1-3 chain reveals distinct differences above as well as below the 
MIT. In retrospect, this interrelation applies also to VgOn because the displacements of the dioxide-like 
1-3-5 chain point along the local z-axis rather than those of the 2-4-6 chain. 

Decomposing the V 3d ti g partial DOS, we observe in figures[2]andEDfor the d x 2 _ y 2 states a similar 
behaviour as in VO2. Due to strong a-type V-V bonding along the vanadium chains the chain center sites 
show distinct splitting of the d x i_ y i DOS. In the case of atom V3 a broad peak is observed from roughly 
— 0.7eV to 0.3 eV and a sharp peak at 1.0 eV. Because of weaker V4-V4 bonding the splitting is smaller 
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Fig. 16 Same representation as used in fieure fT5l but for the vanadium sites WHY A. 



for the V4 site with peaks at approximately — 0.3eV and 0.8 eV. The chain end atoms VI and V2 are less 
affected by in-chain metal-metal bonding and consequently display a more compact d x 2_ y 2 DOS. At low 
temperatures the behaviour of both metal chains changes substantially. Increased V- V dimerization yields 
enhanced bonding-antibonding splitting of the V3 d x 2_ y 2 DOS. Additionally, by virtue of reduced V1-V3 
overlap along c pm t the VI d x 2_ y 2 states tend to localize. Accordingly, the VI d x 2_ y 2 DOS sharpens and 
develops a pronounced peak near the Fermi energy at low temperatures. While the V4-V4 distance grows 
by 0.1 A below the transition, the coupling of the V4 atoms to the chain end sites increases. The reduced 
V2-V4 bond length causes stronger bonding-antibonding splitting of the V2 d x 2 _ y 2 DOS. 

Modifications of the vanadium d xz partial DOS at the transition are less significant. They consist mainly 
of an energetical up and downshift of the states as observed for atoms VI and V2, respectively. Such shifts 
are due to increased and decreased displacements perpendicular to c pru t. Moving vanadium atoms laterally 
off the centers of gravity of their oxygen octahedra modifies the overlap between V 3d and O 2p orbitals. 
The oxygen sublattice itself is hardly affected by the structural transformations. Therefore the shifts of the 
d xz states can be understood in terms of a zigzag-type antiferroelectric distortion of the vanadium chains. 
As the distortion is reduced for the chain center sites V3 and V4, see the former discussion, the effects are 
barely visible in the DOS. In the prototypical case of VO2 lateral displacements are found only in the low 
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temperature monoclinic phase, whereas V4O7 reveals them similar at all temperatures. Since metal-metal 
bonding is less important for d xz orbitals no bonding-antibonding splitting is visible. 

The situation is more complicated for d yz orbitals, which are considerably influenced by different kinds 
of hybridization. Recall from our earlier investigations that these electrons can be involved in metal-metal 
overlap across octahedral faces both along a prut and 6 pru t. Moreover, antiferroelectric-like displacements 
of the vanadium atoms lead to increased overlap between V 3d and O 2p states, which causes energetical 
upshift of the antibonding vanadium dominated states. In the case of V4O7 the atoms V3 and V4 display 
neither relevant antiferroelectric shifts nor V-V interaction parallel to the a pru t-axis. There are simply no 
nearest metal neighbours in this direction. As a consequence, the pronounced double peak structure of the 
V3/V4 d yz DOS indicates bonding-antibonding splitting due to V1-V3 and V2-V4 overlap, respectively, 
which is true at high and low temperatures. Analogous with the finite metal 3-chains of VgOn along o pru t 
we are confronted with vanadium pairs giving rise to VI- V3 and V2-V4 bonding across shared octahedral 
faces. In contrast to the central vanadium atoms the chain end atoms VI and V2, apart from participating 
in V1-V3/V2-V4 overlap along 6 pru t, take part in V1-V2 bonding along a prut . The partial VI and V2 d yz 
densities of states thus consist of features reminiscent of the sites V2/V3 and V1/V4, respectively. Being 
subject to two different types of metal-metal overlap, the d yz orbitals of the chain end atoms undergo two 
intertwining bonding-antibonding splittings, which lead to complicated triangular DOS shapes, especially 
in the low temperture phase, and resemble the shape of the V 3d a\ g DOS of V2O3. In addition, the form 
is affected by the discussed zigzag-type displacements of the chain end atoms, which intensify the metal- 
oxygen bonding. As a consequence, the center of gravity of the V1/V2 d yz DOS shifts to higher energies 
below the MIT. Changes in the shape of the V2 d yz DOS yield a closer resemblance to the VI d yz DOS, 
which can be traced back to an improved coupling between the chain end atoms at low temperatures. As 
the relative displacements of the sites VI and V2 along c prut are reduced, the V1-V2 bond penetrates the 
octahedral faces along a prut more favorably and the overlap becomes more effective. 

In conclusion, the electronic structure of V4O7 is strongly influenced by the coordination of the metal 
sites, as previously reported for VgOn. Generally, while the overlap of V 3d and O 2p orbitals places the 
V 3d t2 g states close to the Fermi energy, the details of the electronic structure are subject to metal-metal 
bonding. In particular, one finds a distinct bonding-antibonding splitting of the d x 2_ y 2 states at the chain 
center sites. In the low temperature phase the same applies to atom V2. All these splittings are analogous 
with the corresponding orbitals in VO2 and VeOn, where they likewise are enforced by V-V dimerization. 
In contrast, the VI d x i_ y i DOS reveals a compact structure and a distinct sharpening below the MIT as 
these states hardly take part in the dimerization along the metal chains. Hence they may be susceptible to 
electronic correlations. For the d xz and d yz states we find noticeable response to the displacements of the 
metal atoms perpendicular to c pru t- As the atoms move relative to the octahedral centers, states with d xz 
and d yz symmetry display energetical up and downshifts in the vanadium and oxygen dominated regions, 
respectively. However, the d yz DOS is affected mainly by the V-V bonding across shared octahedral faces. 
Chain center sites are subject to bonding only along 6 pru t an d show a characteristic two peak structure. In 
contrast, atoms VI and V2 additionally take part in d yz -type overlap along a prut , leading to densities of 
states similar to the a\ g DOS of V2O3. The latter orbital mediates V-V bonding along c nC x = aprut- 

Similar to the VgOn findings, the LDA results indicate small charge transfer between the two inequiva- 
lent metal chains. Recall the limitations of assigning charge to specific atomic sites as discussed in section 
[5] According to table[5]the V2 sites comprise 0.2 electrons less than the VI sites in the high temperature 
modification of V4O7, while V3 and V4 show similar values. Note that in the pairs V1/V2 and V3/V4 the 
sphere radii are almost equal. The interpretation of the low temperature results is a little more complicated. 
However, taking into account the various sphere radii one can state reduced charges for both V2 and V4, 
where the charge deficit is 0.2 electrons per site. Of course, the calculated values of 2.6-2.8 V 3d valence 
electrons do not agree with the ionic picture predicting 1.5 electrons per metal atom - but due to covalent 
V-O bonding this deviation is not surprising. In comparison to VeOn the atomic spheres of V4O7 com- 
prise slightly more charge, which reflects the position of the compound in the Magneli series. We find no 
significant differences in the charges of the particular oxygen spheres since each of them contains about 
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four 2p electrons. Further electrons enter the interstitial empty spheres and thus are difficult to classify in 
vanadium and oxygen-type contributions. Altogether we observe a charge transfer of about 0.1 electrons 
per site from the 2-4 to the 1-3 chain. In the high temperature phase only the chain end sites are involved 
in the transfer. As in the case of V 6 On the different electron count is related to the crystallographical in- 
equivalence of the two vanadium chains. While the 2-4 donator chain becomes slightly more V02-like, the 
1-3 acceptor chain to some extent approaches the V2O3 configuration. These LDA findings qualitatively 
confirm x-ray results by Marezio et al. [10], who applied an ionic picture and reported chains of V 3+ and 
V 4+ ions running parallel to the c pru t-axis. Although an ordered arrangement of varying valence charges 
is confirmed by the calculation, the charge discrepancy is much smaller due to covalent bonding. Nuclear 
magnetic resonance findings by Gossard et al. [9] support chains of differently charged atoms along c pm t 
but indicate incomplete differentiation into 3+ and 4+ valences. The authors observed a slightly increased 
charge ordering below the MIT, which is confirmed by the LDA results. 

As in other Magneli phases, particularly in V 6 On, the electronic structure of V4O7 is a mixture of fea- 
tures characteristic of either vanadium dioxide or sesquioxide. The local atomic environment of the chain 
end site VI bears close resemblance to the V2O3 structure as it is subject to an increased separation from 
its vanadium neighbours at the MIT. We are thus in a position to transfer some of the above conclusions 
for V4O7 to V2O3. First, V-V bonding across octahedral faces along ct pru t seems to have the major effect 
on the shape of the d yz DOS of V4O7 and thus on the shape of the ai g DOS of V2O3. However, also the 
V-V hybridization across the faces parallel to 6 prut is large since it leads to strong splitting of the electronic 
states. Second, due to the reduced metal-metal interaction across the octahedral edges parallel to c pru t the 
vanadium d x 2_ y 2 states are substantially localized and undergo further localization at the transition. 

This supports previous work by Dernier [41], who concluded from a comparison of pure/doped V2O3 
and of Cr203 that the metallic properties of V2O3 are connected to V-V hybridization across shared octa- 
hedral edges rather than to hopping processes in the vanadium pairs along Ch cx - From this point of view 
the MIT arises from structural changes occuring at the transition. Via strong electron-phonon coupling the 
structural distortions translate into a narrowing of the bands perpendicular to the Ch ox -axis. Note that the 
decoupling of the vanadium atoms parallel to the c prut -direction is even stronger for the sesquioxide than 
for V4O7 since the corresponding bond length shortens by 0.12 A instead of 0.06 A. The band narrowing 
eventually paves the way for increased influence of electronic correlations in the insulating configuration. 
So far our treatment of V2O3 concerned only the symmetry breaking PM-AFI transition, which is accom- 
panied by substantial modifications of the interatomic V-V spacings. However, Pfalzer et al. [42] reported 
on local symmetry breaking in the PI phase of aluminium doped V2O3, giving rise to similar structural and 
electronic properties of the AFI and PI phase at least on a local scale. Due to the local distortions one may 
expect a closely related scenario for both the PM-AFI and the PM-PI transition. In general, the reduction 
of the dispersion perpendicular to Ch ox seems to play a key role for the MIT of V2O3. 

7 Systematic aspects of the Magneli series 

Finally, addressing the systematics inherent in the Magneli series we are interested in effects of the chain 
length. Variation of the length systematically modifies the ratio of dioxide and sesquioxide-like regions. In 
this context especially a comparison of V7O13 and V 8 Oi5 is useful because the materials are structurally 
very similar. Nevertheless, the former compound is the only member of the Magneli series which does not 
undergo an MIT as a function of temperature. In a comparison of all Magneli phases, VO2 and V2O3 take 
an exceptional position as they either miss chain end or chain center metal sites. In addition, they display 
higher symmetry than that of the triclinic space group PI {C}). This is likewise true for monoclinic V3O5. 
The latter compound, although even closer related to the sesquioxide than V4O7, is therefore not suitable 
for a detailed understanding of the electronic states in V2O3. Both V5O9 and V 8 Oi5 are characterized by a 
superstructure doubling the unit cell at the MIT. While the metal chains of low temperature V5O9 [74,75] 
lose their inversion symmetry, there are four inequivalent chains for V 8 Oi5. The notation of the vanadium 
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Fig. 17 Notation of the vanadium sites in V7O13 and VsOis, where we observe metal chains of length n = 7/n = 8 
along Cp r ut. In both cases eight inequivalent vanadium sites give rise to two types of chains. The spacial arrangement 
of the chains is similar to that of V4O7 and VeOn, see figures [8] and [5] The table gives measured V-V distances (in 
A) for both the metallic (298 K) and the insulating (20 K) phase of V7O13/V8O15 [60]. Below the transition the latter 
compound develops a superstructure with four inequivalent chains. 



sites in V7O13/V8O15 is given in figure^l F° r the former compound the chains are labeled V1-V3-V5- 
V7-V5-V3-V1 and V2-V4-V6-V8-V6-V4-V2, whereas for V 8 Oi 5 the innermost site is doubled, leading 
to the configurations V1-V3-V5-V7-V7-V5-V3-V1 and V2-V4-V6-V8-V8-V6-V4-V2. Extra sites in the 
insulating superstructure are indicated by tildes. As in V4O7 and VgOn, the oxygen sublattice of VgOis 
hardly changes at the phase transition. The dominating structural modifications are again concerned with 
the V-V dimerization along c prut , see the bond lengths in figure [171 

The study of V7O13/V8O15 is again based on LDA calculations, which were performed analogous with 
the previously discussed Magneli phases. The partial DOS shows three groups of bands corresponding to 
O 2p, V 3d t2g, and V 3d e a g states. Due to a reduced formal electron count of 1.29 and 1.25 electrons per 
metal site for the n — 7 and n — 8 compound, respectively, the filling of the t2 g group is slightly smaller 
than for V4O7 and VeOn. Site-projected partial densities of states for V7O13 and VgOi5 show the typical 
features discovered before. In particular, the shape of the DOS curves is completely understood in terms 
of the vanadium coordination, which is true for all Magneli compounds. Hence it is useful to concentrate 
on representative atomic sites: the central site V7 and the last but one site V3 of the 1-3-5-7 chain. 

Figure^]shows site-projected partial densities of states for the sites V3 and V7 arising from the room 
temperature structures of V7O13 and VgOi5. The representation of the V 3d symmetry components refers 
to the local rotated reference frame with the z and x-axis parallel to the apical axis of the local octahedron 
and the c pmt -axis, respectively. V-O overlap combined with crystal field splitting results in two groups of 
antibonding states. While the t 2g dominated region reaches from about — 0.6eV to 1.8 eV, one observes 
additional t 2g contributions above 2.0 eV, where the e a g states dominate. In consistence with the positions 
of V7O13 and VgOi5 within the Magneli series, the latter are fairly small, which resembles the behaviour 
of VO2 rather than of V2O3. Contributions due to O 2p orbitals in the energy interval of figure^] which 
arise from covalent V-0 bonding, comprise less than 10% of the total DOS at the Fermi level; they are not 
included in the figure. As a matter of fact, the in-chain V-V distances in the high temperature modifications 
of V7O13 and VgOis are very similar, see figurefTTl Of course, additional V7-V7 and V8-V8 bonds occur 
for chains of length n = 8, which are not available for chains of length n = 7. The oxygen sublattice of the 



32 



U. Schwingenschlogl, V. Eyert: The vanadium Magneli phases V n 02 n -i 



DOS (1/eV) 
3 



V 7 13 at 298 K 



V3 site 




V7 site 




1 2 

(E - Ep) (eV) 



1 2 

(E - E F ) (eV) 



DOS (1/eV) 



V 8 15 at 298 K 



V3 site 




V7 site 




1 

(E - E F ) (eV) 



1 

(E-Ep)(eV) 



Fig. 18 Site-projected partial V 3d t2 g densities of states (DOS) per vanadium atom for high temperature V7O13 and 
V8O15. The orbitals refer to the local rotated reference frame. Sites V3 and V7 reveal little discrepancy. 



compounds is almost identical, which likewise applies to the antiferroelectric-like lateral displacements of 
the metal atoms perpendicular to c pmt and the rotation of the metal chains off the c plut -axis. Compared to 
V4O7 and VgOn the latter effects are weaker due to the dioxide-like nature of the long chain compounds. 
The importance of the local atomic coordination for the DOS shape suggests a similar electronic structure 
for both metallic V7O13 and VgOi5, which is confirmed by figure[T8l 

The local d x 2_ y 2, d yz , and d xz densities of states resemble one another even with regard to the details of 
their shapes, which applies not only to the sites V3 and V7 but also to all other metal atoms. Due to strong 
er-type metal-metal bonding along c pmt the chain center atoms are characterized by a distinct splitting of 
the d x 2_ y 2 DOS. In the case of V7 we hence find sharp peaks at about — 0.3 eV and 1.1 eV. Because of the 
increasing V-V distances at the chain ends the bonding-antibonding splitting is noticeably smaller for the 
V3 d x 2_ y 2 DOS, where peaks occur at about 0.2 eV and 0.7 eV. As a consequence of less efficient 7r-type 
V-V overlap via the d xz orbitals the respective DOS is rather compact. In contrast, a substantial splitting 
in bonding and antibonding branches can be found in the d yz DOS. For sites V3 and V7 this shape traces 
back to metal-metal interaction along 6 piu t- The strength of the splitting is roughly identical for both sites. 
To conclude, the site-projected densities of states of high temperature V7O13/V8O15 are well understood 
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V 8 15 at 20 K 




(E - E F ) (eV) (E - E F ) (eV) 

Fig. 19 Site-projected partial V 3d t2 g densities of states (DOS) per metal atom for sites V3 and V3 of low temperature 
VgOis. At high temperatures these sites are equivalent. The orbitals refer to the local rotated reference frame. 



from the discussions in the preceding sections. By virtue of very closely related metal-metal bond lengths 
in the compounds similar DOS curves arise. Thus it is surprising that only one of the materials undergoes 
an MIT. While the charges in the 1-3-5-7 chain of V7O13 and the 2-4-6-8 chain of VgOis are more or less 
equally distributed among the metal sites, the other chains show a small shift of charge towards the chain 
centers, but generally such charge transfer effects are very small. In the calculation the metallic phases of 
V7O13 and VgOis are barely distinguishable as far as the local electronic behaviour is concerned. 

The investigation of the low temperature VsOi5 modification reveals the same mechanisms found pre- 
viously in the Magneli series. Depending on the bond lengths to the adjacent vanadium sites along c pru t 
one finds either increased splitting or enhanced localization of the d x 2_ y 2 orbitals. The V-V bond lengths 
indicate strong V-V dimerization in the 1-3-5-7, 2-4-6-8, and 2-4-6-8 chain, but not in the 1-3-5-7 chain. 
While the innermost V-V distances of both the 2-4-6-8 and 2-4-6-8 chain elongate at the MIT from 2.87 A 
to 2.98 A and 3. 14 A, they shorten for the other chains from 2.78 A to 2.61 A and 2.57 A. Thus the dimer- 
ization patterns are inversed. As a consequence, in the low temperature structure all vanadium sites, except 
for VI, V3, V5, and VI, are involved in strong V-V overlap along c plut . The increased interaction leads to 
a stronger bonding-antibonding splitting of the d x 2_ y 2 DOS. For atoms affected by V-V dimerization one 
additionally obtains energetical upshifts of the d yz and d xz densities of states. In accordance with earlier 
findings they show typical changes in the electronic structure known from the MIT of VO2. Although the 
energetical separation between the d x 2_ y 2 states and the remaining e y states of VgOis grows in the low 
temperature phase it is still far from complete. The described Peierls-like behaviour does not account for 
sites VI, V3, V5, and VI, which are characterized by decreasing V-V overlap below the MIT. Due to the 
reduced metal-metal bonding the d x 2_ y 2 densities of states sharpen and develop a pronounced peak close 
to the Fermi level. Because of their sesquioxide-like coordination the behaviour is not surprising for sites 
VI and VI. However, the separation from neighbouring vanadium atoms forces V3 and V5 to display the 
same localization of the d x 2_ y 2 orbitals. Hence electronic correlations may be as important for these states 
as for the V203-type chain end sites. The question why one of the four vanadium chains in the low temper- 
ature VsOi5 superstructure prefers localization to dimerization is difficult to answer. The calculated V 3d 
valence charges do not give rise to any systematics. Partial £23 densities of states are depicted in figure[T9l 
For the dioxide-like V3 site we recognize both an increased splitting of the d x 2_ y 2 DOS and an energetical 
upshift of the d yz and d xz curves. In contrast, the structurally separated V3 site is characterized by a very 
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compact d x 2_ y 2 DOS located almost completely above the Fermi level. In comparison to the sesquioxide- 
like VI site of V4O7, see figure[T5l the filling of the low temperature d x 2_ y 2 states is significantly smaller 
due to the reduced electron count. Summarizing, these findings support our MIT scenario as a mixture of 
features typical of either vanadium dioxide or sesquioxide, which we developed in the previous sections. 

We are nevertheless left with the question why the compound V7O13 does not undergo an MIT. From 
purely structural considerations the main difference between V7O13 and VsOis is the amount of dioxide- 
like center sites. The local environments and thus the electronic properties of the individual sites are very 
similar. However, we point out that the MITs in the Magneli series arise as combinations of an embedded 
Peierls instabilty inherent in the metal chains and electronic correlations relevant for the chain end atoms. 
Since compounds with small vanadium chains are dominated by chain end sites and thus more susceptible 
to correlations, their MITs are most likely driven by sesquioxide-like states. On the other hand, the Peierls 
mechanism should be dominating for materials with long chains. This scenario of different driving forces 
of the MITs is consistent with the behaviour of the transition temperatures in the Magneli series as seen in 
figure^ With only one exception (n = 5) the temperatures decrease monotonously from V3O5 to V7O13, 
which does not show a transition. In the case of V5O9 the peculiar superstructure evolving below the MIT 
might cause the observed deviation. Transferred to our picture this kind of systematics corresponds to the 
decreasing influence of chain end sites as compared to the growing impact of the chain centers. When the 
chain length increases the influence of the local electronic correlations at the chain ends diminishes. At the 
same time the Peierls-like mechanism becomes stronger. With an increasing number of chain center atoms 
available for V-V dimerization the energy gain resulting from such a modification becomes more and more 
advantageous. With the experimental transition temperatures in mind we thus attribute the MIT of V§Oi5 
to a Peierls-like scheme; for V7O13 neither of the two mechanisms is strong enough to induce a transition. 

8 Conclusions 

The vanadium Magneli phases V„02 n -i form a homologous series of compounds with crystal structures 
comprising dioxide and sesquioxide-like regions. Studying the MITs in this material class thus paves the 
way for a more complete understanding of the phase transitions of the prototypical oxides VO2 (n — ► 00) 
and V2O3 (n — 2). Analyzing those V 3c? <2g states which are involved in the phase transitions allows us 
to discuss the changing influence of electronic correlations and electron-lattice interaction on going from 
VO2 to V2O3. We have developed a unifying description of the various crystal structures including both 
the rutile structure of VO2 and the corundum structure of V2O3. Our systematic representation is based 
on a regular three-dimensional network of oxygen octahedra partially filled by metal atoms, giving rise to 
characteristic metal chains of length n. Due to the comprehensive picture of the structures, the electronic 
bands can be grouped into states behaving either VO2 or V203-like. The detailed electronic structure and 
consequently the MITs of the Magneli phases turn out to be strongly influenced by the local metal-metal 
coordination. The phase transitions arise as results of both electron lattice interaction within the dioxide- 
like and electronic correlations within the sesquioxide-like regions of the crystal. Dioxide-like vanadium 
sites reveal the characteristic features of the embedded Peierls instability responsible for the MIT of VO2. 
The combination of dimerization and antiferroelectric-like displacements of the metal atoms by means of 
strong electron-lattice interaction causes splitting of the d x 2_ y 2 = dn states and energetical upshift of the 
dyzldxz states. Sites related to V2O3 are characterized by strongly reduced V-V overlap. Thus the d x 2_ y 2 
states become more localized. Consequently, correlation effects play a more important role for these states. 
In conclusion, the electronic structures and the metal-insulator transitions of the Magneli phases remain a 
crucial test case for theories aiming at a correct description of both VO2 and V2O3. 
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